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ABSTRACT 


This report culminates the work accomplished during a three year 
design project on the automation of an Environmental Control and 
Life Support System (ECLSS) suitable for space travel and 
colonization. The system would provide a comfortable living 
environment in space that is fully functional with limited human 
supervision. A completely automated ECLSS would increase astronaut 
productivity while contributing to their safety and comfort. The 
first section of this report. Section 1.0, briefly explains the 
project, its goals, and the scheduling used by the team in meeting 
these goals. Section 2.0 presents an in-depth look at each of the 
component subsystems. Each subsection describes the mathematical 
modeling and computer simulation used to represent that portion of 
the system. The individual models have been integrated into a 
complete computer simulation of the C0 2 removal process. In 
Section 3.0, the two simulation control schemes are described. The 
classical control approach uses traditional methods to control the 
mechanical equipment. The expert control system uses fuzzy logic 
and artificial intelligence to control the system. By integrating 
the two control systems with the mathematical computer simulation, 
the effectiveness of the two schemes can be compared. The results 
are then used as proof of copcept in considering new control 
schemes for the entire ECLSS . Section 4.0 covers the results and 
trends observed when the model was subjected to different test 
situations. These results provide insight into the operating 
procedures of the model and the different control schemes. The 
appendix, section 5.0, contains summaries of lectures presented 
during the past year, homework assignments, and the completed 
source code used for the computer simulation and control system. 
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1.0 INTRODUCTION 


1.1 Project Description 


For prolonged missions into space and colonization outside the 
earth's atmosphere, development of Environmental Control and Life 
Support Systems (ECLSS) are essential to provide astronauts with 
habitable environments. ECLS systems for Space Station Freedom 
(SSF) require semi-autonomous operation to allow environmental 
control without constant supervision by crew members. The Kansas 
State University Advanced Design Team is in the process of 
researching and designing a control system for an ECLSS like that 
on Space Station Freedom. 

The ECLS system for Freedom is composed of six subsystems. The 
Temperature and Humidity Control (THC) subsystem maintains the 
cabin temperature and humidity at comfortable levels. The 
Atmosphere Control and Supply (ACS) subsystem insures proper cabin 
pressure and partial pressures of oxygen and nitrogen. To protect 
the space station from fire damage, the Fire Detection and 
Suppression (FDS) subsystem provides fire sensing alarms and 
extinguishers. The Waste Management (WM) subsystem compacts solid 
wastes for return to earth, and collects urine for water recovery. 
Carbon dioxide and other dangerous contaminants are removed from 
the air by the Atmosphere Revitalization (AR) subsystem. The Water 
Recovery and Management (WRM) subsystem collects and filters 
condensate from the cabin to replenish potable water supplies, and 
processes urine and other waste waters to replenish hygiene water 
supplies . 

At this time, automation and control of these subsystems have not 
been fully developed or integrated. A fully integrated and 
automated ECLS system would increase an astronaut's scientific and 
observational productivity as well as contribute to their safety 
and comfort. 


1.2 Three Phase Design Schedule 


Kansas State University implemented a three phase approach to 
facilitate the design of a control scheme for an ECLS system. Each 
phase, consisting of one academic year, represented an evolution 
and advancement of previous progress. 

The first phase consisted of information gathering and determining 
the particular tasks required for design of the ECLS system. This 
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accumulated knowledge led to the present organizational structure 
centered on six interconnected subsystems. 

The second phase examined the Air Revitalization subsystem. The 
concept of a series of mathematical models providing input to a 
control system was chosen. Prototype models of the C0 2 Removal 

Assembly governed by a crude expert system controller were 
developed. 


The third phase concentrated on refining the C0 2 Removal Assembly 
and comparing two control schemes. The two control systems 
compared are a classical proportional-integral-differential 
controller and an expert system fuzzy logic controller. The 
purpose of this study is to enhance the knowledge of these control 
approaches so choices can be made for the control scheme for the 
entire ECLS system. 


1.3 Third Year Goals 


the proposed goals for the third phase of the design 
proj ect were to combine the control systems of the six subsystems 
and form an overall control system for ECLSS with fault diagnosis. 
However, due to lack of control systems for the individual 
subsystems, the goals for phase three were reevaluated. 

The overall objective of the final year is to develop and compare 
expert and classical systems of control on a computer simulation of 
the C0 2 Removal Assembly of the EC LS S . 

Goals for reaching the final objective begin with creating a 
mathematical model and a computer simulation of the C0 2 Removal 
Assembly. Concurrently, development of the classical and expert 
systems of control were performed. The next goal is to integrate 
the control systems and the computer simulation together and 
evaluate and compare the effectiveness of each control system. The 
comparison will be used as a proof of concept to evaluate the use 
of expert systems to control the entire ECLSS. 

A list of the goals for the third and final year are as follows: 


1. Complete the computer simulation of the C0 2 Removal Assembly. 

2. Create a set of rules for the expert control system of the 
C0 2 Removal Assembly. 

3. Create a classical controls system for the C0 2 Removal 
Assembly. 
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4. Establish a means of communication between the mathematical 
model and the two controls systems. 

5. Analyze the dynamic response of the simulation and compare the 
two methods of control. 


1.4 Academic Year Time Table 


The year started with an introduction to the advanced design teams 
objectives for the project. Several lectures given by faculty and 
graduate members of Kansas State University introduced the design 
team to mathematical modeling, simulation and control. This 
introduction lasted until September 25 . 

The next step was to plan objectives for the first semester, and 
decide what should be .accomplished for the third year. A 
comparison of expert system controls and classical system controls 
for a subassembly of the ECLSS was decided upon as the third year 
main objective. 

Three modeling groups and two controls groups were formed to 
develop models for the individual parts of the C0 2 Removal 
Assembly. The modeling of the desiccant beds, the blower/ 
precooler, and the C0 2 sorbent beds began about October 23 r , with 
completion deadlines planned for November 10 . These models were 
to be integrated together forming a computer simulation of the 
overall process. A presentation of the progress was given on 
November 2 5 th . 

Documentation of the semesters work started on December 2 n and a 
semester report was submitted to the faculty advisors on January 

_ * rd 

23 . 

The final semester goals were to refine the math models formulated 
during the fall semester, complete implementation of controls on 
the C0 2 Removal System, and create a user interface using X 
Windows . 

On January 30 th two modeling groups were formed. One to refine the 
math model of the sorbent beds and another to find information on 
the inputs and outputs to the C0 2 Removal Assembly . Two other 
groups were formed to implement classical and expert controls on 
the C0 2 Removal Assembly. On March 4 , work on a Graphical User 
Interface (GUI) using X Windows was begun. All phases were 
completed for the All University Open House on April 4 with 
displays in the College Engineering and the College on Arts & 
Sciences. 
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Documentation and a presentation for the USRA/NASA Design Program 
was begun on April 6 fc * and continued until the day of the conference 
on June 15 th . 


1.5 Design Team Description 


The Advanced Design Team at Kansas State University is composed of 
students from several academic disciplines. Currently 
participating disciplines include computer science, and mechanical 
engineering and chemical engineering. The team's Graduate Teaching 
Assistant is an electrical engineer. Plans are under way to 
recruit students from electrical and computer engineering for the 
final semester. Faculty support comes from the mechanical, 
electrical, chemical, and computer engineering departments as well 
as the computer science department. 
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2.0 MATHEMATICAL MODELING 

2.1 C0 2 Removal Assembly 

2.1.1 Introduction 


The Carbon Dioxide Removal Assembly, designed to remove carbon 
dioxide from the cabin air, involves removal of C0 2 by molecular 
sieves. The process is required to remove carbon dioxide generated 
by th f J es P lrat ory processes of the astronauts and to maintain 
acceptable levels of carbon dioxide within the cabin. 

Figure 2 . 1-1 is a block diagram representation of the C0 2 Removal 
Assembly. The system takes input air from the Temperature Humidity 
Control Subsystem (1), and valves (2,11) direct the air flow 
allowing it to flow across one of the desiccant beds (3,10) which 
dehum idify the air using Zeolite 13X and Silica Gel. The moisture 
must be removed to avoid poisoning the desiccant found in the 
adsorbing sorbent bed (8,14) . Because the dry air is heated in the 
process, it rs forced across a heat exchanger (6) by a blower (5) 
anc * air is cooled before being sent through a sorbent bed. The 
sorbent beds remove the carbon dioxide by means of Zeolite 5A 
which acts as a molecular sieve adsorbing the carbon dioxide. The 
dry air returning from the molecular sieves through unidirectional 
control valves (13,9) is revitalized by the moist desiccant of the 
second desiccant bed (10). After the air is rehydrated it is then 
returned to the Temperature and Humidity Control Subsystem (12) and 
redistributed throughout the cabin. 



Figure 2.1—1: co 2 Removal Assembly 


5 









Concurrently, a second desorbing sorbent bed (14) is being heated 
causing the separation of the carbon dioxide from the desiccan . 
The desorbed carbon dioxide is drawn from the bed by means of a 
pump (16) and is sent to an accumulator tank (18) . After the 
adsorbing desiccants have become saturated, the desorbing beds are 
once again dry. The control valves (5,7,15) redirect air flow in 
the system. The previously adsorbing beds begin the desorbing 
process and the previously desorbing beds begin adsorbing. The 
system is presently configured to cycle every thirty minutes. 

Mathematical models of the various components were created to allow 
analysis of the subassembly's performance. The role of the 
modeling is to duplicate the actual systems response to a given set 
of parameters. Knowing how an actual system should respond, it is 
possible to explore control systems for use in governing the 
subassembly. The control systems regulate the state variables 
throughout the subassembly. 


2.1.2 Desiccant Beds 


DESCRIPTION 

The purpose of the desiccant beds is to remove water vapor from the 
incoming air. This function is necessary because water will poison 
the zeolite used in the C0 2 adsorption process . Water vapor 
removal is achieved by means of two desiccants, Silica gel and 
Zeolite 13X. High-humidity air coming into the C0 2 Removal 
Assembly flows first over the Silica gel, which removes some of the 
water and brings the relative humidity down to a low lf ve l* Since 
Zeolite 13X works well at low relative humidities, it is then used 
to remove most of the remaining water vapor before the air is sent 
on to the blower and precooler. The exiting air is not only dry, 
but also heated from the release of energy required to condense the 
water vapor. As adsorption is taking place in one bed, the other 
desiccant bed is rehumidifying the air returned from the C0 2 
sorbent beds through a desorption process. The desorbing cycle is 
iust the reverse of adsorbing in that hot, dry air flows across the 
Zeolite 13X first, then over the Silica gel, and is cooled and 
humidified in the process. 


MATH MODEL 


Assumptions 

1 Equilibrium relative humidity is a linear function of load for 
Silica gel, while that for Zeolite 13X can be approximated 
with two straight lines. 
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2 . 


All heat transfer is between the air and water only. 

3. Water in the desiccant is evenly distributed at all times. 

4. Beds remain at room temperature. 

5. Air is heated before desorption. 

6. The water is removed and released at constant pressure. 

7. The specific heat of the air is constant. 

8. Enthalpy of condensation and vapor saturation pressure are 
accurately represented by linear and fourth order least 
squares curve fits, respectively. 

9. Equilibrium relative humidity, which is a function only of the 
load on the desiccant bed, is achieved. 

Equations 

A thermodynamic analysis of the air flowing through the bed yields 
the following equations used for the mathematical model. Equation 
1 is a curve fit for vapor saturation pressure as a function of 
saturation temperature. Values for the curve fit were obtained 
from the steam tables in an appendix of Thermodynamics ,, — An 
Engineering Approach by Cengel and Boles. The expression is 


p sat = ■ 3972 + . 0629T+ . 001099T 2 +1 .105xl0~ s T 2 +6 . 192x10 1 T* . (1) 

From the same text, relationships equating pressures, mass, and 
relative and absolute humidities are shown in equations 2 , 3 , and 
4 as 

P v=4> P sat' 

_ .622 P v 


( 2 ) 

(3) 


The law of mass conservation can be applied to the model when 
evaluating the air inside the desiccant bed. The mass of dry air 
in the bed (m a ) equals the total mass of air in the bed (mf) minus 
the mass of water vapor in the air (hO state as 

m a =m f -m v . (S) 
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The load on the desiccant can be defined as the mass of water vapor 
absorbed versus the mass of desiccant in the bed expressed as 


L= 


m 


ads 


m 


tank 


( 6 ) 


Equilibrium load curves were provided by Dr. Byron Jones of Kansas 
State University from an ASHRAE reference. Data points taken from 
these curves of the load versus relative humidity for the 
desiccants were curve fit using a least squares method fortran 
program written by Dr. Kirby Chapman, professor of Mechanical 
Engineering at Kansas State University. The resulting linear curve 
fits are given in Equations 7, 8, and 9. The fit for Zeolite 13 X 
was approximated using the two straight lines of equations 8 and 9. 
The results are 

<b= — - — (Silica gel) , (7) 

v .5263 


<J>= . 4 L (Li£. 17, Zeolite 13X) , (®) 


$=. 068+40 (L-. 17) (L> . 17 , Zeolite 13 X) . (9) 

The mass of water removed from the air can be determined using the 
thermodynamic relationship 


m z =m v -(nm a • 


( 10 ) 


Enthalpy, the incoming air temperature, and the outgoing air 
temperature are related by the thermodynamic reations 

h fff =2S02~2 . 389 T in , UD 


and 


T = T • + 
A 2 M in 


( m f -m z ) Cp ' 


( 12 ) 


The remaining equations are simply relationships for the rates of 
mass absorbed or desorbed, and the change in mass of air vapor and 
air in the bed versus time. These equations govern the mass 
transfer of the water from the desiccant and the air. Utilized in 
a finite time step process, the following expressions determine the 
success of the desiccants in removing water from the air. 
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( 13 ) 



dt dt 


dm 


dt 



dm f _ dm z 
dt dt 


( 14 ) 


( 15 ) 


The symbols used in the mathematical model are defined as follows. 

p sat = vapor saturation pressure (kPa) 

T = temperature of air being evaluated (°C) 

P v = partial pressure of the water vapor (kPa) 

<p = relative humidity 

u = absolute humidity 

p = total pressure of the air (kPa) 

mv = mass of vapor in the air (kg) 

m f = total mass of air in the bed (kg) 

m a = mass of dry air in the bed (kg) 

L = load on the desiccant 

m ada = mass of water adsorbed by the desiccant (kg) 

"“tank = mass of desiccant material (kg) 

m r = mass of water removed from the air (kg) 

h £g = enthalpy of condensation (kJ/kg) 

T in = temperature of incoming air (K) 

T 2 = temperature of outgoing air (K) 

c p = specific heat of air (kJ/kg*K) 

dm.,! . = rate at which the desiccant adsorbs water (kg/s) 

dt 

dm r = rate at which water is removed from the air (kg/s) 
dt 

dm- = change in mass of vapor in the air (kg) 
dt 

dnu = change in total mass of air in the bed (kg) 
dt 


MODELING TECHNIQUES 

A computer program was written to simulate the performance of the 
desiccant beds in time. This was accomplished by choosing a small 
time step, and then evaluating the above equations for the mass in 
the bed during that time step. The first calculations on this mass 
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determine the composition of the air based on the input conditions. 
Next, the equilibrium relative humidity is found from the load on 
the first desiccant to come in contact with the air flow — Silica 
gel for adsorbing, Zeolite 13X if desorbing. The change in the 
amount of water in the air is found from the difference between the 
input and equilibrium states, with removal of water from the air 
considered positive. The temperature change of the air is 
dependent on the amount of water removed. However, changing the 
temperature of the air also alters its relative humidity (but not 
the absolute humidity) so that it no longer matches equilibrium. 
Therefore, an iterative procedure is required to find the point 
where the output temperature and its corresponding humidities are 
consistent with the equilibrium relative humidity and the new 
amount of water in the air. Execution of the program showed that 
an average of five iterations were needed. With this done, the 
mass of air is sent on to the other desiccant, and the calculations 
are repeated to produce the final output conditions of the air from 
the bed. 

It should be noted that little distinction is made between the 
adsorbing and desorbing cycles. This is because the desiccant 
itself is unaware of the intended function; it merely reaches 
equilibrium with the conditions it is given. Naturally, adsorption 
will occur when cool, humid air passes over desiccant with a low 
loading. In order for the loaded desiccant to be desorbed by the 
dry air on the return trip, the air must first be heated by means 
of a heat exchanger because hotter air holds more water vapor. 
Since the current system does not account for heating the air, the 
program sets the input temperature to a sufficiently high value on 
the desorption cycle. 


RESULTS 

Many simulations were run with the program, each time varying one 
parameter to observe its effect. Since the temperature and 
relative humidity of the air being sent to the C0 2 desorption 
process (and then back to the cabin) are the most important, it is 
no surprise that they have the greatest effect on the performance 
of the desiccant beds. Temperature is the most influential of all 
parameters because relative humidity is a function of temperature. 
The mass of the desiccant is also an important parameter, since it 
affects the loading. Ideally, any problems with the amount of 
absorbed water vapor could be solved by varying the desiccant mass 
within the bed. 

The following are expected normal operating conditions, inputs, and 
experimentally determined values: 

Mass flow rate of air .2 kg/s 

Adsorption input temperature 300 K 
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Desorption input temperature 363 K 

Adsorption input relative humidity 100% 

Total pressure of the air 101.325 kPa 

Specific heat of air 1.006 kJ/kg*K 

Initial mass of water adsorbed • 0 kg 

Mass of each desiccant 40 kg. 

Figures 2.1-2 and 2.1-3 show the behavior of the desiccant beds 
from startup under standard conditions. One initial condition of 
the process assumes the desiccants to be completely dry. 
Desiccants have a tendency to retain some moisture which will not 
be retrievable during the desorption cycle. Over time, the amount 
of residual water will asymptotically approach an equilibrium value 
which will represent the maximum amount of irretrievable water 
vapor retained by the desiccant. This characteristic accounts 



Time (hours) 


Figure 2.1-2: standard Temperatures Leaving Desiccant Bed 
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0.3 



Figure 2.1-3: Standard Humidities Leaving Desiccant Bed 

for the transient responses of the adsorbing and desorbing cycles 
in both figures. Notice from the figures that after startup, the 
output temperatures and humidities reach a steady pattern. After 
reaching a steady state, the average outputs of the adsorption 
cycle tends to be a temperature of about 63 °C with roughly 6% 
relative humidity, while the desorption cycle is returning air at 
averages of about 58 °C and 20% relative humidity. 

Because heat is generated when water vapor is adsorbed by 
desiccant, the amount of water vapor adsorbed is directly related 
to the output temperatures from the beds. In addition, the amount 
of water vapor adsorbed by a desiccant during a thirty minute cycle 
may be represented by a decaying exponential function. Examining 
the first 30 minute adsorbing cycle of Figure 2.1-2, one may notice 
a falling output temperature due to the decrease in heat generated 
by the decaying amount of water vapor adsorbed during each time 
step. 

Since the adsorption and desorption cycles are essentially 
reciprocal, the exponential decay of water vapor released during 
the desorption cycle, and the increase in output temperature can 
also be explained. The desorption cycles are conducted with the 
air flow being at a higher temperature. That fact coupled with the 
removal of C0 2 result in the differing magnitudes of change in 
relative humidity. 
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2.1.3 Blower and Precooler 


DESCRIPTION 

The blower/precooler is the second process of the C0 2 Removal 
Assembly. This process utilizes a variable speed blower to force 
cabin air through the C0 2 Removal Assembly and a crossflow heat 
exchanger to cool air received from the desiccant beds. Air 
leaving the precooler is then directed on to the sorbent beds where 
carbon dioxide is removed. 


HATH MODEL 

Assumptions 

1. Pressure drop across the cooler is negligible. 

2. Water specific heat and air density are constant properties. 
Specification of Heat Exchanger 

1. Heat exchanger effectiveness is 0.80. 

2. Heat exchanger coolant is water. 

3. Mass flow of coolant is 3.79 kg/min (500 lbm/hr.). 

4. Inlet temperature of coolant is 15 °C (59°F) . 

5. Cross-sectional area is 11.1 m 2 (120 FT 2 ). 

Equations 

Equation 1 is the maximum heat transfer rate that can be drawn from 
the air by the heat exchanger. Equation 2 is the actual heat 
transfer rate using the heat exchanger effectiveness. The concepts 
give 


2 ™ 


(t)=c . 

lax ' — 7 — p, ail 


M 


air 


(— h, i ( — ) 


( 1 ) 


^(t^eQ^t). 
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Equation 3 is the temperature of the air leaving the precooler as 
a function of time given the inlet temperature, the actual heat 
transfer rate and the mass and specific heat of air. Symbolically 
this can be written as 


T 

— h , out 


(t) =T. 


■h,i 


(t) -- 




C . M 

— p, air — air 


( 3 ) 


The inlet temperature function used to test the program is given as 


T u . (t) =2t +T . 

— h , i ' — 9 — — D 


( 4 ) 


This function varies with time because the outlet temperature from 
the desiccant beds is not constant. Note .that the equation used 
does not represent the actual desiccant bed outlet temperature, but 
is only an approximate linear function used to test the response of 
the precooler model at extreme conditions. 


The symbols used in the above equations are defined as follows: 


c 

v -'p, air 
6 

M air 

Qact(t) 

Qmax(t) 

T c .i 

T hfi (t) 

Th , out (^) 

T 0 

t 


specific heat of air (kJ/kg*K) , 
heat exchanger effectiveness, 
mass flow rate of air into cooler (kg/s) , 
actual heat transfer between two fluids (kW) , 
maximum amount of heat transfer between two fluids, 
(hot inlet air and cold coolant water) (kW) , 
inlet temperature of cooling water (K) , 
inlet temperature of hot air (K) , 
outlet temperature of cooled air (K) , 
initial temperature of inlet air (300 K) , 
time (minutes) . 


RESULTS 

Figure 2.1-4 compares the inlet and outlet temperatures of the 
precooler. As shown in the figure, the inlet temperature of the 
air increased from 300 K to 420 K over 60 minutes while the 
precooler outlet temperature varied between 290 K and 313 K. This 
precooler temperature range was used for sorbent bed inlet air 
temperatures in other modeling for that component. 
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Temperature In Temperature Out 


Figure 2.1-4: Precooler Inlet and Outlet Temperatures 


2.1.4 C0 2 Sorbent Beds 


DESCRIPTION 

The C0 2 Sorbent Beds are the third step of the C0 2 removal process. 
Their purpose is to separate C0 2 from the air, return the air to 
the desiccant beds and send the removed C0 2 gas to an accumulator 
tank. This is done using two adsorption beds containing Zeolite 
13X. 


MATH MODEL 

Assumptions 
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1. The exiting air temperature equals the temperature of C0 2 in 
the bed. 

2. The heat transfer coefficient of Zeolite 5A is a linear 
function of temperature. 

3. Power supplied to the beds is either on (1000 J/s) , off (0 
J/s) , or being removed at 1000 J/s. 

4. Thermal equilibrium for the sorbent beds negates dependence on 
bed length. 

5. Assume that a four man loading supplies C0 2 at an average rate 

of 5.046x 10* kg/s (this can be a function of time). 

Equations 

The mass of C0 2 contained within an adsorbing bed at a time T is 
equal to the mass adsorbed at some time T-t plus the mass 
transferred during time t. The result is of the form given in 
Equation 1 below. This mass transfer is brought about by the 
sorbent bed removing the C0 2 based on a difference in the 
equilibrium and actual partial pressures. This mass transfer is 
given by equation 2 shown below. The expressions are 

m T =iri T _ t +m t: *At, (1) 


and 




( ^OOj ^ equilibrium 


) V , 


tank 


R 


CO , ± tank 


( 2 ) 


The ideal gas law applied in the above equations provide a simple 
relationship between the mass flow rates that are desired and the 
partial pressures which are known. 

The energy involved in the mass transfer and accompanying phase 
change results in the bed temperature being raised. In addition, 
during the desorbing phase the bed is heated to drive off the C0 2 
and this results in a further increase in bed temperature. This 
temperature change is governed by the following energy balance: 


( 3 ) 


Finally, the equilibrium partial pressure curve was derived from 
curve fitting data provided from Marshall Space Flight Center. 
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This data shows that the equilibrium partial pressure of the bed 
and air flow is a function of both the corresponding temperatures 
and most importantly the bed loading, or current level of absorbed 
mass. This information is expressed as 

PPC0 2 =. 1333 exp (10 2oad -3.065 + . 02622 3^-4 . 684E-5 *I’ bed 2 ) . ( 4 ) 


The symbols used in the above equations are defined as follows: 

Tbed = temperature of sorbent bed (K) , 
mcoz = mass of C0 2 transferred (kg) , 
p c02 = partial pressure of C0 2 (kPa) , 

Pequii = equilibrium C0 2 concentration of bed (kPa) , 

RC 02 = C0 2 gas constant (kPa m 3 )^(kg K) , 

Vbed = volume of sorbent bed (m ) , 

HEAT = power applied to bed (J/s) , 

m air = mass of air in bed (kg) , 

m bed = mass of sorbent bed (kg) , 

mats = mass of C0 2 absorbed in bed, 

cv air = specific heat of air (J/kg K) , 

cv C 02 = specific heat of C0 2 (J/kg K) , 

cVb e d = specific heat of sorbent bed (J/kg K) . 


RESULTS 

The system is designed to include two sorbent beds which alternate 
between the adsorbing and desorbing roles. While one bed is 
adsorbing the C0 2 flowing through the system, the other is being 
heated and its previously adsorbed C0 2 is released and pumped out 
to the accumulator tank. The Figure 2 . 1-5 shows the loading curves 
for the two beds. The increasing curve is indicative of the bed 
that is loading, while the decreasing bed's loading is shown as the 
curve that is falling off. 

The work done on the subroutine involved a total overhaul of the 
previous semester's model due to unacceptable limitations . in the 
earlier version's performance. This work included enhancing the 
accuracy of the model's portrayal of the actual phenomena, and 
increasing the subroutine's compatibility with the main program. 

After the fundamental flaws were corrected, the problem of fine 
tuning the desorption process was examined. Two major criteria 
were established as defining the problem. First it was neccessary 
to desorb all the C0 2 in the half hour cycle, and second, it was 
neccessary to provide almost pure C0 2 gas to the accumulator tank 
which feeds the Bosch reactor. 
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SORBENT BED 

Bed Loading (KG/KG) 



Figure 2.1-5 C0 2 Bed Loads 


The first problem was solved by incorporating a heating/cooling 
jacket to the sorbent bed. This allowed the temperature of the bed 
to be raised which resulted in a lower affinity for the adsorbed 
C0 2 . The immediate problem with this was the need to cool the bed 
before returning it to the adsorbing cycle. A 1000 watt 
heating/cooling jacket was found to be adequate to accomplish both 
of these ends. 

The result of this bed heating was that the C0 2 gas was desorbed 
into the bed to mix with the residual air still in the tank after 
the cycle switch. During the first three minutes of the adsorbing 
cycle the bed is vented into the Temperature and Humidity Control 
assembly to avoid contamination of the input C0 2 gas for the Bosch 
Reactor. The final result of the heating and cooling curves can be 
seen in Figure 2.1-6. 

The adsorbing process is not unidirectional in that the bed 
achieves an equilibrium not neccessarily in phase with the desired 
result. The temperature dependency of the equilibrium partial 
pressure results in a transient desorption phase in the beginning 
of the adsorption cycle. The resulting change in equilibrium 
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SORBENT BED 

Bed Temperature 



Figure 2.1-6 CO z Bed Temperatures 


nar1 . ial nressure results in the mass transfer of C0 2 as the bed 

Figure S 2 t^e 

the bed desorbs and the pump removes C0 2 . As the bed heat P th^ 

g^ph 1 b ^owe v e r ^"a s *t he^ l^d* *de sen: bs* i t s & CO 2 'its^equi librium partial 

pressure tails off, • /i"f %£&?* £2E,£?^i2™ 
bed temperature and the th at as the bed returns to the 

further. The end result . . , nressure is very low and 

adsorbing cycle, its equilibrium partial J>n Y 

it is able to immediately begin adsorbing C0 2 . 
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1.4 


SORBENT BED 

Bed Partial Pressure (kPa) 



Figure 2.1-7 CO z Pressure for Adsorbtion Cycle 


2.1.5 Pump and Accumulator 


DESCRIPTION 

The C0 2 pump and C0 2 accumulator tank is the f ina re leased 

CO, Removal 2JS Vhf ac^ulato^TaS 

stores^the * the «* * edUCti °" 

Assembly. 


MATH MODEL 

Assumptions 

1. The C0 2 pump is a fixed displacement, rotary vane pump 

2. Pump is 100% efficient. 
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3 . 


pvunp operates under adiabatic and isentropio conditions 

4. Accumulator tank is perfectly insulated. 

5. co 2 is an ideal gas. 

6. C0 2 specific heat is constant. 


Equations 

The rate of CO, mass pumped to the 

determined by the following equa fso=rev/sec. ) and the volume 
displaced “p^re^luS ™ « » obtain the rate 

equation as 


ih= 


P_sp_V 
R T 


( 1 ) 


The temperature on the outlet s 1 ®* sid e P 'of’ th^pump* and the 
^um P p r tirrhat d ^”prTssu“e Ch to S tS pump is constant. This 
gives 


T = T .(—£) (1 - 1/w • 
1 po ■ L pi ' p . 


( 2 ) 


The time change in 
tank is given by 


enthalpy at both the inlet and outlet of the 

H in =m p ((T po -27 3)C p +h 0 ) , < 3 > 


and 

H out ='»»<< I 'r- 273 )C n +h «> • 


( 4 ) 


Sele^'knTwing^e 

volume. The relationships are given by 

m Ti = V T P Ti / (RT Ti ) , (5) 


and 


U Ti =(T Ti ~213)C v m Ti . 


( 6 ) 
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The overall mass within the tank as a function of time is found by 
subtracting the C0 2 drawn from the tank from the mass pumped into 
the tank and adding it to initial C0 2 . This can be stated as 

m T m o ^ cp * • ^ ^ ^ 


The current tank temperature and pressure is found using the first 
law of thermodynamics and the ideal gas law respectively. These 
concepts give 


T T ~T T +ni C o i ^C Vi co 1 '^i/ ( ^co 2 , i +m T^ ' 


( 8 ) 


and 


P rp Tp/ Vp • 


(9) 


The symbols used in the mathematical model are defined as follows 


m p 

sp 

V P 

Ppi 

R 

Tpo 

Tpi 

P T 

k 

H ln 

H 0 

Cp 

H ou t 

m 0 

m Ti 

V T 

P T i 

T Ti 

U Ti 

m T 

fstp 

U T 

T t 

Cv 

P T 


= mass flow rate of C0 2 from the pump (kg/sec) , 

= speed of the pump (revolutions/sec), 3 

= volume displaced by the pump per revolution (m ) , 

= inlet pressure to pump (kPa) , 

= ideal gas constant (kPa*m /kg*K) , 

= outlet temperature of the pump (K) , 

= inlet temperature of the pump (K) , 

= pressure of C0 2 in the tank (kPa) , 

= specific heat ratio of C0 2 , 

= enthalpy of inlet C0 2 stream to the accumulator (J/s) , 
= enthalpy of C0 2 at the reference temperature (J/kg) , 

= constant pressure heat capacity of C0 2 gas (J/kg*K) f 
= enthalpy of outlet C0 2 stream from accumulator (J/s) , 
= mass flow rate of C0 2 leaving accumulator (kg/sec) , 

= initial mass of C0 2 in accumulator tank (kg) , 

= volume of accumulator tank (m ) , 

= initial pressure inside tank (kPa) , 

= initial temperature in tank (K) , 

= initial internal energy of tank (J) , 

= current mass of C0 2 in the accumulator tank (kg) t 
= time elapsed between calculations (sec ) , 

= current internal energy of tank (J) , 

= current temperature of tank (K) , 

= constant volume heat capacity of C0 2 (J/kg*K) , 

= current pressure inside tank (kPa) . 


RESULTS 

The development of a pump and accumulator simulation was assigned 


22 


as homework to the design team as an introduction to math modeling. 
The homework summary and source code can be found in the appendix 
under section 5.4 Homework 1 Summary as well as results for the 
mathematical model . The graph of conditions is in Figure 2 . 1-8 



Temperature Pressure 


Figure 2.1-8: Homework 1 Case 1 Data 

The linear curve corresponds to tank pressure, and is listed 
against the right hand axis. The upper curve is the tank pressure 
and is read off the left axis. The results show tank pressure 
approaching a maximum value while the pressure continues to 
increase. 


2.2 C0 2 Reduction Assembly 


INTRODUCTION 

The following description is used to define the input variables to 
the C0 2 removal model. The C0 2 Reduction Assembly consists of a 
Bosch reactor, heater, compressor, condensing heat exchanger and 
dynamic separator, and recuperative heat exchanger. A schematic of 
the C0 2 Reduction Assembly is given in figure 2.2-1. 


DESCRIPTION 
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Figure 2.2-1: C0 2 Reduction Assembly 


The C0 2 flows from the C0 2 accumul^ ^ c0 is mixed with an 
control valve to the Bosch Rea * Assembly. The two 

excess of Hydrogen from X ^hich is collected in the Carbon 

react to give solid carbon, , whi J reac tion is less than 10 % 
Filtration Cartridge, and wat ® • recycled. After passing 

efficient so the remaining separator, the 

dr ied^ r e t ant^ a^reTm i xed^w iJoming £2 

&£* S£ Z 99% P of S the CO^ i s ^reduced to ’solid Ca rt on 

and water. 

The reaction that takes place within the reactor is: 

CO. + 2H-, = 2H 7 0 + C + HEAT * 


suppUed d by Y Z ZZZzZV S? thT^ct^n^nevIr 
generates enough heat to be eelf-sustainmg. 


CO z REDUCTION DESIGN 

1. Bosch maximum capacity for 1 unit- 4 astronauts or 
C0 2 per day. 


8.8 lb of 
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2 . 


Modes of operation of the C0 2 Assembly 


a. Normal mode - The Bosch is reducing C0 2 to water & 
Carbon (S) . The single pass efficiency is less than 10%, so 
the reactants have to be recycled and passed back through the 
Bosch. Less than one percent of the exiting reactants are not 
reduced to water & solid Carbon. The Bosch will operate for 
90 days before servicing the Carbon Filtration Cartridge is 
necessary. 

b. Standby mode - Everything is powered & ready to go except 
all valves are closed and the compressor is off. 

c. Shutdown mode - The heater & compressor are off and all 
valves are closed. All sensors are working. 

d. Purge mode- The system is being purged with nitrogen. The 
purge is drawn off through the nitrogen purge/bleed vent. The 
compressor & heater are off. 

e* Unpowered mode - No electrical power is applied to system. 

3 . Process startup- The process starts in the unpowered mode and 
switches to the shutdown mode. The system is checked for 
leaks and if none are detected the Carbon Dioxide Reduction 
Assembly is purged with Nitrogen. While the Assembly is being 
purged, the heaters in the Bosch are turned on and kept at a 
constant 200 F for two hours. This is to drive off any 
moisture accumulated in the Bosch during servicing of the 
Carbon Filtration Cartridge. After the two hours the leak 
check and purge are finished and the heater temperature is 
increased to 1050 °F. The compressor is started and the purge 
Nitrogen is circulated around the system. When the reactor 
temperature reaches 1050 °F the reactants are introduced. The 
average time from leak check to normal operating mode is 12 
hrs . 
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Metabolic C0 2 

Desian Point 

Ranqe 

Flow Rate,lb^day 
Temperature, F 

8.8 

8.80-17.60 

70 

60-85 

Pressure , psia 

18 

14.7-20 

H 2 Feed 



Flow Rate, lb^ day 
Temperature, F 

.80 

. 80-1.60 

75 

75-100 

Pressure , psia 

30 

14.7-30 

Product Water 



Flow Rate,lb^day 
Temperature, F 

7.20 

7.20-14.40 

60 

60-90 

Pressure , psia 

30 

14.7-30 

Bleed 



Flow Rate,lb^day 
Temperature, F 

1.12 

1.12 

75 

65-90 

Pressure , psia 

18 

14.7-20 

Electric Power 



28 VDC, W 

341 

306-606 

115 AC, W 

186 

170-3120 

Heat Rejection,W 



To Air 

529 

494-818 

To Coolant 

238 

181-461 


Table 2.2-1 Flow Rates of Inputs 2 

1 "Carbon Dioxide Reduction Description", Boeing Aerospace & 
Electronics. Huntsville, Alabama. April 20, 1990, Doc # 2-H8RG-RJK- 
198. 

2 "Bosch Carbon Dioxide Reduction (Bosch - III) Subsystem", Volume 
1, Life Systems Inc. October 1987. 
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2.3 Temperature and Humidity Control 

Subsystem 


INTRODUCTION 

Because the C0 2 Removal Assembly model takes air from the 
Temperature and Humidity Control (THC) Subsystem, it was necessary 
to research the THC Subsystem and determine the effects it may have 
on the air entering the C0 2 Removal Assembly. 


DESCRIPTION 

The principal function of the Temperature and Humidity Control 
(THC) Subsystem is to maintain a comfortable environment for the 
astronauts in Space Station Freedom (SSF) . Air in the cabin is 
continuously circulated through the THC system at 340 cubic feet 
per minute. Temperature and Humidity are controlled by a 
condensing heat exchanger and a slurper, respectively. Depending 
on the temperature change needed to keep the cabin air within 
certain specifications, a temperature control valve determines the 
amount of air passed through the condensing heat exchanger. Air 
not passed through the heat exchanger is bypassed to the exit side 
of the heat exchanger. As air passes through the heat exchanger, 
water vapor is condensed and drawn into a slurper to remove excess 
water vapor in the air and prevent excessively high cabin humidity. 
After leaving the heat exchanger, the air is. pulled across a mixed 
flow fan, and passed through an air straightener before it is 
returned to the cabin. 


EFFECTS ON C0 2 REMOVAL ASSEMBLY 

According to available THC documentation, air drawn by the C0 2 
Removal Assembly is taken from the THC system immediately after the 
temperature control valve, and returned just before it reaches the 
air straightener. This configuration is not acceptable for the 
following reasons. First of all, the air drawn by the C0 2 removal 
assembly cannot be taken after the temperature control valve 
because in some cases, the control valve may bypass all air flow 
around the heat exchanger leaving none available to the C0 2 Removal 
Assembly. Secondly, because the air returned from the C0 2 Removal 
Assembly may not be within the established cabin parameters, it 
should not be returned after the THC air conditioning process. 

The model will take air from the high pressure side of the THC fan 
and return it to the THC Assembly at the inlet. Two assumptions 
will correspond to this configuration. Because air is taken from 
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the THC after the heat exchanger, inlet conditions into the C0 2 
Removal Assembly will be assumed to remain within cabin parameters. 
Secondly, the flow rate of air from the CO z Removal Assembly will 
be assumed much smaller than that being drawn from the cabin so the 
retruned air cannot cause a false response in the temperature 
control valve. 


CABIN PARAMETERS 

The cabin pressure operates between 14.5 and 14.9 psia, and the 
temperature will be kept between 65 and 80 °F. Relative humidity 
and partial pressure of C0 2 will be maintain within 25 to 7 5 % 
humidity and 3 to 12 mm of mercury, respectively. 

The only cabin parameter which will act independently of the THC 
system is the production of C0 2 . A schedule which approximates the 
production of C0 2 by the astronauts for a 24 hour period is 
described in Section 2.4. 

The effects of these parameters on the model will be tested within 
and outside of the ranges given. 

2.4 Cabin Model 


DESCRIPTION 

A model of the cabin was produced to simulate the effects of C0 2 
production and removal on the cabin atmosphere. The model 
simulates temperature, pressure, and relative humidity levels 
within the cabin by three different functions: a constant value, 
sinusoidal and step functions varying within specified parameters. 
For each time step, the model evaluates the amount of C0 2 produced 
within and removed from the cabin and determines the current 
partial pressure of C0 2 inside the cabin. 


MATH MODEL 

The model allows for the cabin conditions of temperature, pressure, 
and relative humidity to be simulated in several ways. Relative 
humidity and temperature can be varied by use of either a sine or 
step function and will fluctuate between any given parameters 
establish within the program. The model also allows for varying 
the cycle time of each function. Because normal cabin pressure 
conditions only fluctuate between 99.9 and 102.7 kPa, the pressure 
is only simulated by either a constant value or a sine function. 
After the values for the cabin pressure, temperature, and relative 
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humidity are determined for a particular time step, the program 
evaluates the amount of C0 2 produced by the astronauts for that 
time step. 

Figure 2.4-1 shows the C0 2 production within the cabin per 24 hour 
period. Beginning at midnight, time zero, the C0 2 production 
levels are at the lowest value because the astronauts are sleeping. 
At 0600 hours, all four crew members awake and for the next four 
hours each takes an hour of exercise, producing the highest level 
of C0 2 procudtion. From 1000 to 1400 hours, normal breathing 
processes are maintained as the crew performs the daily duties, and 
from 1400 to 1800 hours each memeber again takes an hour of 
exercise. At 1800 hours normal activities are resumed until the 
astronauts go to bed at 2200. 
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Figure 2.4-1: Cabin C0 2 Production 


Next, the model calculates the mass of C0 2 in the cabin by taking 
the previous amount of C0 2 in the cabin, adds the mass produced and 
subtracts the amount removed for that time step. The mass of the 
air at the current cabin conditions is determined using the ideal 
gas law in Equation 1 and stated as 

m _ (CPress) ( CVol ) , 1} 
aiz (R) (CTemp ) 


The symbols used in the above equations are defined as follows: 

M air = mass of air in the cabin, 

CPress = current cabin pressure. 
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CVol = cabin volume (101.115 m 3 ) , 

R = gas constant (8.314/29), 

CTemp = current cabin temperature ('C). 


The mole numbers for the cabin air and the C0 2 in the cabin are 
calcuated by dividing the mass of the air and C0 2 by the molecular 
mass of air and C0 2/ respectively. Next, the mole fraction of CO? 
to air is calculated by dividing the number of C0 2 moles by the 

Sf r, m °iS S ° f air, 4 - and . the P aritial pressure is calculated by 
multiplying the current cabin pressure by the mole fraction. The 

caculated cabin partial pressure of C0 2 is then checked to 

determine if the C0 2 removal assembly needs to be turned on to 
remove any excess C0 2 . 
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3.0 CONTROLS 


3.1 Classical Control 


DESCRIPTION 

The C0 2 removal sub-assembly is responsible for maintaining the 
partial pressure of C0 2 with in normal limits as the astronauts and 
other equipment and experiments produce it. NASA grades air 
quality by the partial pressure of C0 2 , with normal C0 2 pressure 
being 0.0667 kPa. When the C0 2 partial pressure is above 0.4 kPa 
the air is classified as "degraded" and above 1.015 kPa the 
condition is classified as "emergency". The C0 2 removal sub- 
assembly removes C0 2 from the cabin environment and stores it as a 
gas in a C0 2 accumulator tank until the Bosch reactor breaks it 
down to solid carbon and water. 

The C0 2 removal sub-assembly uses a variable speed fan to force air 
through the system 1 s beds , ducts and heat-exchangers . The 
desiccant beds and the C0 2 sorbent beds operate on 30 minute 
cycles, where one bed adsorbs mass for 30 minutes while the 
companion bed is desorbing. After 30 minutes the beds reverse 
roles and the full adsorbing bed desorbs its mass while the empty 
desorbing bed adsorbs mass. 

CLASSICAL CONTROLS 

There are two inputs that control the operation of the C0 2 removal 
sub-assembly, the partial pressure of C0 2 in the cabin and the 
pressure of C0 2 in the C0 2 accumulator tank. The cabin C0 2 pressure 
input is used as input to a classical control to maintain the cabin 
C0 2 pressure. If the partial pressure of C0 2 in the cabin deviates 
from the desired 0.0667 kPa the system would modify the air flow 
rate. 

The input from the C0 2 accumulator tank was based on the gas 
pressure in the tank. The Bosch reactor is an important producer 
of fresh water and a shortage of C0 2 may mean a corresponding 
shortage of fresh water. The Bosch reactor shuts down if the 
pressure of the supply C0 2 (the C0 2 tank) dips below 101.125 kPa, 
so the system is turned on if the pressure in the C0 2 accumulator 
tank drops below 137 kPa. This safety buffer of 3 6 kPa assures 
that the tank pressure should not go below the lower limit of 
101.125 kPa. 

Internal to the C0 2 removal sub-assembly are controls that maintain 
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the pressure of the C0 2 accumulator tank and a valve that is 
positioned before the C0 2 accumulator tank and after the C0 2 pump 
that controls the purity of the C0 2 entering the tank. 


The cabin air is driven through the system by a variable speed, 
zero-inertia fan that is controlled to maintain cabin pressure of 
0.0667 kPa. Classical control of the fan speed is accomplished by 
using a proportional-integral-differential (PID) compensator in a 
negative feedback loop. The PID compensator uses an error function 
6, defined as the difference between the actual C0 2 cabin pressure 
and the desired cabin pressure. The magnitude of the change in the 
pump speed is given as 


a. . = 5 + -^+/Sdt. 

^fanspeed 


(i) 


The fan speed is then adjusted by this amount, increasing or 
decreasing the tank pressure. 

The valve between the C0 2 accumulator tank and the C0 2 pump serves 
two purposes. One is to direct C0 2 gas to the accumulator tank when 
the pressure in the desorbing C0 2 sorbent bed is within 1% of the 
equilibrium pressure of C0 2 for the bed. This insures that the gas 
that is directed to the C0 2 tank is almost entirely C0 2 . The other 
purpose is to quickly evacuate the air from the C0 2 sorbent bed 
that just switched to the desorbing cycle. As the beds switch, the 
full bed that is just beginning to desorb contains cabin air and 
C0 2 trapped in the absorbent material. For the first several 
minutes of the desorbing cycle the gas removed from the bed is air 
and, as the pressure in the bed decreases as the air is removed, 
the temperature of the bed increases and the equilibrium pressure 
of the C0 2 trapped in the Zeolite begins to increase. As the 
pressure of the bed and the C0 2 equilibrium pressure converge, the 
purity of the C0 2 gas leaving the bed increases. When there is a 
difference of greater than 1% between the pressures, the valve 
directs the gas back to the exit gas from the adsorbing C0 2 sorbent 
bed and turns the C0 2 pump to its maximum speed to expedite the 
emptying of the bed. 
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3.2 Expert Systems Control 


IMPLEMENTATION 

The simulation of the Carbon Dioxide Removal Assembly can be 
controlled by an expert system written in CLIPS using fuzzy logic. 
The simulation for the physical system is written in FORTRAN. The 
purpose of using FORTRAN is that an existing FORTRAN simulation has 
already been developed by mechanical engineers of the NASA group. 
Last semester we struggled with choosing between C and FORTRAN as 
a simulation language. The simulation equations were taken from 
the existing FORTRAN simulation and implemented in C. The C model 
then communicated with CLIPS to make a separate model apart from 
the classically controlled simulation. 

Originally the computer science students felt that it would be 
easier to integrate CLIPS into the C environment. Since they would 
be implementing the CLIPS program into the simulation, they felt 
they should work with a program that was most familiar to them, 
hence C. Later on when the expert control was running, it was 
found that it was a painfully slow working with the simulation. 
This is caused by too much file I/O overhead. The reason for this 
is that CLIPS cannot communicate or link to any programming 
language other than itself. The problem is sharing variables 
between two languages. On the one hand, we did not want to 
implement the whole model within CLIPS. It is not that easy to 
program a simulation using an expert systems programming language. 
Rules do not get fired in the order that one expects. On the other 
hand, we did not want to implement the whole system in C either. 
Programming a recursive expert systems controller in C can be quite 
a struggle. It would be easier to use an expert systems program 
that was designed to do just that. Therefore, we were left with 
the job of integrating the two into one environment. Initially we 
used a file sequencer that monitored the reading and writing of the 
variable files between C and CLIPS. This was extremely slow and 
used about 90% of the processor of a Soulbourne SPARC computer 
running UNIX System Release V. 

Since the original design was slow and the system administrators of 
the computer resources weren't happy that it took so much processor 
time, we ventured to design a new system. Semaphores was one of 
the solutions that were brought up, but no success on 
implementation was ever achieved. A semaphore is a process that 
gets "forked” off from the initial process. What a semaphore does 
is protect what is known as a critical section. In our case, the 
critical section is the file being passed between CLIPS and C. 
This file contains all of the variables pertinent to the running of 
the simulation. Some of the variables simply get read by the 
processes, and other get changed in the process. However, most, if 
not all, variables in this file get changed at one time or another. 
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The problem is that not both processes can be writing to this file 
at the same time. This would cause chaos. So this is where a 
semaphore becomes useful. A semaphore would allow the C program to 
start its simulation process, and when it was ready, "fork" off the 
expert control CLIPS process. You then tell the semaphore where 
the critical section of code exists (the reading and writing of the 
variables file) , and then both processes can run simultaneously 
waiting for the next one to hand control to it. Again, a semaphore 
in this case would keep the variable file from being written to at 
the same time by two different processes. However, as stated 
before, a working implementation of semaphores was never worked 
out. 

So then we move onto the third and final design which we are 
currently using. The C simulation was dropped and the FORTRAN 
simulation was used for both the expert control and the PID control 
to maintain a completely consistent environment. 

At this point we are in the Spring 1992 semester of implementation 
of the project. We have decided to use FORTRAN as the simulation, 
C and Xwindows for a user interface, CLIPS for the expert systems 
controller, and FORTRAN for the PID controller. We have also 
decided to use a Soulbourne SPARC UNIX workstation as our platform 
of choice. Our reasons behind this are simple; it is an extremely 
fast computer, it multitasks, *and hard disk space on this system is 
plentiful . 

The engineers programmed the FORTRAN simulation and PID controls 
while three computer scientists programmed the expert systems 
control and the Xwindows user interface. Much of the semester was 
spent by the engineers getting a working "bug free" simulation 
running so the controllers could be employed. The expert systems 
controller was built within the first month and then modules were 
stubbed for tests. After this, efforts were placed on getting the 
user interface to work. 

There is not much to the expert systems in terms of lines of code. 
However, this does not mean that it does not do a lot. Some of the 
expert system was designed like the PID controller because there 
was really no expertise that could be used in making a decision. 
Situations such as a valve require only two positions, i.e. on or 
off. Other such devices require only If/Then statements, and no 
fuzzy logic was used in determining what variables to change. For 
the devices that can take advantage of fuzzy logic, CLIPS becomes 
very powerful. Once the engine to determine membership functions 
has been written, it can be used over and over to control many 
different devices. All that need be added to it are one or two 
lines at the beginning of the CLIPS code that describe what ranges 
the variables should exist in. The engine takes care of the rest. 
One such example might be as follows; 

(fuzzy tank-pressure low -150 300 400), 
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(fuzzy tank-pressure high 300 400 650), 

(fuzzy cabin-pressure low -100 .058 .075), 

(fuzzy cabin-pressure high .058 .075 100). 

Presently, the only component within the C0 2 Removal Assembly which 
is controllable is a variable speed fan used to draw mass through 
the sorbent beds, which remove carbon dioxide from the atmosphere. 
The controller monitors the Cabin C0 2 partial pressure to determine 
when the pump speed should be adjusted to maintain safe C0 2 levels. 


PRESSURE CONTROL BY EXPERT SYSTEM METHODS 

The expert system uses triangular functions to control the 
simulation. A triangular function consists of three values: low, 
medium and high, as shown in Figure 3.2-1. 


Fuzzy Logic 



Figure 3.2-1 Fuzzy Logic Membership Triangle 


A function is used to calculate a percentage belief when the value 
being considered is in the range low to high. When the value does 
not lie in the range low to high, the percentage belief is zero. 
A belief is calculated with Equation 1 when the value being 
considered is in the range low to mid, and Equation 2 is used when 
the value is in the range mid to high. Given by 
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value -low 


( 1 ) 


; / 
mid -low 


and 


high-value 
' ♦ 
high -mid 


( 2 ) 


The percentage belief is used to directly determine the amount of 
change that must be made. This expert system uses two triangles to 
control the simulation. The left triangle represents the low 
pressure function. The right triangle represents the high pressure 
function. There is also overlap between the high and low 
triangles. This is not uncommon in fuzzy logic. The intersection 
point of the two triangles is chosen to correspond to the target 
control value and to a 50% belief in both triangles. This is done 
so that when the system variable deviates from the target value, 
the belief is immediately greater than 50% in one of the triangles 
which prompts the system to try to correct it. The slope of both 
triangles is adjusted to control the rate at which the expert 
system changes the simulation. Pump speed, pump duration, and 
pressure deviation are factors used in determining the adjustments 
to the triangular functions. The pressure can be controlled more 
accurately when the pump speed is changed more often. However, 
this can cause wear on a pump and must be taken into consideration. 
The definition of the functions in CLIPS are as follows. 

(deffacts start 

(state open) 

(fuzzy temp low 0 258 338) 

(fuzzy temp high 268 348 600)) 

When the percentage belief in a low pressure is greater than 50%, 
Equation 3 is used to adjust the pump speed. Likewise, when the 
percentage belief of a high pressure is greater than 50%, Equation 
4 is used to change the pump speed. In this way the expert system 
is able to control pump speed by monitoring the tank pressure as 
given by 


NewPumpSpeed=01dPumpSpeedx(l+%beliefcold ) , ( 3 ) 


NewPumpSpeed^OldPumpSpeedx (l-%beliefhot) . ( 4 ) 


After trial runs were executed using these equations, it was 
decided to adopt a more fluid control equation. It employs a 
normalized belief, and is less prone to overshoot and repeated 
searching for the desired value. Equation 5 shows the method of 


36 



employing this normalized belief as . 


NewPumpSpeed=01dPumpSpeedx (1+ (2x%belief cold-1 ) ) . (5) 


This improved control equation was then adopted into all full 
simulation exercises. 

The expert system and simulation communicate by writing a temporary 
file. The two programs work in lock step. In other words, one 
program runs one cycle, then the other program runs one cycle. The 
FORTRAN simulation was modified to run only one time step and then 
shell out to the operating system to call CLIPS. The simulation 
will have to read in variables from a file each time that it runs. 
It must then save the variables to the same file after each run. 
The expert system will be acting in the same way with one 
exception. Instead of running continuously, the expert system will 
run only once, make the necessary changes to the file variables, 
and then exit; thus handing control back over the simulation 
program . 
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4.0 DYNAMIC SYSTEM SIMULATION 
4.1 Introduction 


description op the co 2 REMOVAL SUB-AHSEMELY 


The objective of the co , 

fr0nl cabin "ir and store is to remove carbon 

further processing. The sub-a^Ui^i accumulator tank for 

components: two desiccant beSJ ^M Consists of nine main 

sorbent beds, a co 2 pump a co l„ a bl ° wer ' a pre-cooler, two 
schemes. The function of eanh ^ umulator tank, and control 
described below. sub-assembly component is 

cabin air. Thi^is^on^simJltfne Ca ? in alr ' and humidify fresh 
working together at a slt one^i 7 by two desiccant bear 

removing water from us*f cabi"ai? 9 Whil ® one bed is 

water to carbon dioxide free, fresh calin ai?. b6d iS rel eesing 

ss.°sa« 

carbon dioxide from dry cabinet t' a-I^ 1 ° ne bed is removing 9 
stored carbon dioxide to the CO, accumVato? t^nk.^ releasi "9 

bed 6 a'nd^se^on^to^the^ 0 ^ Cabb °" di °~ d ® **« one sorbent 
accumulator tank stores bba accumulator tank. The co 2 bt 

-s. and passes tt&STS. % ~ 

Pf°™atIon S^Jh^cor^ve? in^thTcal, ° f the blower based on the 
methods are used. The first dicrr- Two diff erent control 

and"? a PID approach. The second inlth a S a ola ssical method 
and fussy logic to accomplish lontrlf^er^hew^r SySta " 

HPWr 


development op the simulation program 


a T n h ! control V the^main 3°£“ b ° accurately node! 

? et operating time using a Fortran 2 R ® duction Sub-Assembly 
ulation program is composed of a naiTsSu^^ogrSTand 
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several subroutines that model the main sub-assembly components. 

Before a main program was written, individual subroutines were 
written to model each main sub-assembly component. After the 
subroutines worked successfully on their own, they were 
integrated together using the main program. 

The function of the main program is to initialize all variables 
that are common to two or more subroutines, and pass these 
variables through CALL statements. By initializing most of the 
variables in the main program, any changes to these variables can 
be done by accessing the main program only, and not each 
individual subroutine. Any variables that are common to only one 
subroutine were kept localized within that subroutine. The main 
program operates on a set operating cycle at a constant time 
step. Currently, the operating cycle of the main program is set 
at 24 hours or 84600 seconds, and the time step is set at one- 

tenth of a minute or 6 seconds. The name of the main proqram is 
COOL. FOR 

The function of the subroutines is to calculate common subroutine 
variable values for a given time step. By keeping all 
computations within the subroutines, it is easier to locate, 
assess and adjust erroneous data. The names of the subroutines 
are as follows: 


desiccant beds: DESSBED. FOR 

blower/pre-cooler: BLOWCOOL.FOR 

sorbent beds: SORBED. FOR 

sorbent beds/C0 2 pump: SORPUMP.FOR 

accumulator tank: C02TANK. FOR 


Analysis of the entire C0 2 Removal Sub-Assembly reveals 4 state 
variables for each main component. Given as 

1) mass flow rate, 

2) pressure, 

3) temperature, 

4) relative humidity. 


These state variables are coded and localized for each subroutine 
ky the type of variable it is (mass, pressure, temperature, 
relative humidity) , which subroutine it is in (DESSBED, BLOWCOOL 
SORBED, SORPUMP, C02TANK) , and whether it is at the inlet or 
outiet (in, out). For example in the DESSBED subroutine, mass 
flow m is represented by M for mass, DB for desiccant bed, and 
IN for inlet, yielding MDBIN . Similarly, inlet pressure and 
temperature are PDBIN and TDBIN and outlet mass, pressure, and 
temperature are MDBOUT , PDBOUT, and TDBOUT . Following this 
pattern, the following variables were coded: 
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blower/precooler: MBPIN, PBPIN, TBPIN, MBPOUT, PBPOUT, TBPOUT 

sorbent beds: MSBIN, PSBIN, TSBIN, MSBOUT, PSBOUT, TSBOUT 

pump/tank: MPTIN, PPTIN, TPTIN, MPTOUT, PPTOUT, TPTOUT 

When passing variables from one subroutine to the next, the 
outlet variable for one subroutine will become the inlet variable 
f° r the following subroutine. As a result, mass out from the 
desiccant beds becomes the mass in to the blower/precooler . 
Consequently, the main program passes the previous outlet 
variables in the CALL statements, but receives them with the 
inlet variables in the subroutines. 


4.2 Classical Control Results 


INTRODUCTION 

The simulation with controls needed to be thoroughly tested. 

This would result in two benefits. First it would be possible to 
^ e hermine if the physics of the C0 2 removal process were being 
correctly modelled. Second, it would allow an insight into the 
abilities of both the system and the controllers to handle 
various situations. 

The method used of evaluating the control systems was to 
determine which "weighting factor" provided the most desired 
response. The major characteristic looked for in the solution 
was the ability of the controller to dampen out initial 
transients, and settle upon a closely bound mass flow rate and 
hhsrefore C0 2 rate. This resulted in the system being run at a 
nearly constant rate which greatly reduces wear on the fan due to 
cycling. 


Although many tests were run, the test conditions used for the 
evaluation of the controllers was a simple twin step function 
with an initial offset. It was desired to maintain cabin C0 2 at 
0.0667 kPa throughout the test. The initial value in the cabin 
was set at 0.07 kPa. The C0 2 production rate was initially given 
as 1.7*10 kg/sec, indicative of resting astronauts. At four 
hours into the simulation this value was increased to 7.0*10 5 
^ 9 / se c, a number representing a double sized crew performing hard 
Finally at eight hours the level was decreased to 3.0*10 5 

^ 9 /sec a level appropriate for the standard 4 man crew performing 
typical functions. 

DESCRIPTION 

The classic, or PID, controller was designed around the 
corrective algorithm given by 
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— error 
• • , error dt 
m=m+ ( + 

kl k2 


errordt 


where m refers to the mass flow rate through the blower. In its 
initial form the values of kl, k2, and k3 Sere all equli to 
unity. This resulted m two major effects. First the controller 
was able to quickly vary the flow rate resulting in the 
controller exhibiting a very high frequency. Second the 
derivative s terms influence was very small. The Figure 4 2-1 
shows this controller's response to the test conditionSdetailed 
n the preceding paragraph. The partial pressure of CO, in the 
cabin corresponds to the top curve and is scaled along the riaht 
hand axis. The mass flow rate through the system is the bottom 
curve, and is scaled along the left hand axis. bottom 


CLASSIC RESPONSE 

CREW EXCHANGE (l) 



Figure 4.2-1 System Response with Weighting ( 1 , 1 , 1 ) 

thil fnitfJl ha f 1 i ttle if an y dampening evident, and so 

tnis initial set of constants scored poorly on the scale of 

£• ThiS led t0 the need to increase the impact of the 
rivative term, and also to lower the frequency of the 

controller as the original constants lead to value searching at 
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unrealistic rates. 

a se !r° nd tr ^ tl } e value of kl and k3 were increased to 10. 

This would result m a slower frequency due to the controller 

Svstl^fo^h maS f £- OW at a slower rate ' an <J a better dampened 
system as the relative impact of the derivative term would be 

increased. The results of this controller when subjected to a 
simiiar test are shown in Figure 4.2-2. This controller was able 
. K^ h u SVe ^appreciable amount of dampening during the four to 
eight hour interval corresponding to the highest C0 2 production 

J a * e ‘ ?£T eVer ; a i otl ?er times it was unable to achieve dampening 
and so this set of weighing factors did not represent a 9 

satisfactory solution. 


CLASSIC RESPONSE 

CREW EXCHANGE (id) 
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Figure 4.2-2 System Response with Weighting (.i,i,.i) 


The next attempt was with the value of k2 still at 1 and the 

an i k3 set at 30 • The s y stem response is shown in 
Figure 4. 2 -3 Here we see some significant dampening, and the 
mess flow rate stays relatively well bounded. This controller 
could be labelled as acceptable, but it was decided to see if a 
further improvement could be found. 




TME (HOURS) 


Figure 4.2-3 System Response with Weighting (.03,1,. 03) 


trend would indicate that increasing kl and k3 results in a 
better controller. In pursuit of that trend the next controller 
was run with the weighing factors even larger. The constant k2 
was again left at 1 to provide dampening, while kl and k3 were - 
increased to 50 to reduce the controller's frequency. The 
results from this test were largely similar to the results from 
the controller run previously with a few minor differences, and 
are given in Figure 4.2-4. First the transient spikes in partial 
pressure from the step changes were a little larger, though still 
acceptable. Second this controller though not as capable 
at dampening during the big C0 2 production period, it was a more 
effective controller during the final four hour period. 

There is no reason that the value of kl and k3 had to be left 
equal to each other. Since the system was well behaved and 
smooth, it was not necessary to incorporate a large integral 
term. This fact allows us to assign a very large value to k3 and 
in essence reduce the PID controller to a nearly PD controller. 

By reducing the input from the integral term, it was possible to 
increase the contribution of one of the remaining terms and 
maintain a similar controller. 



CLASSIC RESPONSE 

CREW EXCHANGE (50) 



Figure 4.2-4 System Response with Weighting (.02,l,.02) 


lue i° f k l y as alread y fairly small, it was decided to 
f? « the value of kl back to 25 to increase the effectiveness 
£ho th Proportional term. The net result was a controller with 

do at kl = 25 ' k2 = 1, k3 = 100. These constants 

d° ^°t represent a calculated attempt at optimizing the 

a 1 °g ical qualitative approach to examine the 
effect of the different error terms on the overall responses to 

Figu£l S 4.*2-5 he ^ ltS response to the test case is shown in 

controller exhibits several characteristics. First it 
?hf f o™5 r °? jK lar ? e spike ! n Partial pressure corresponding to 

0 084 kPa Of rn e £? functions. The maximum value attained was 

0.084 kPa of C0 2 . The duration of the spike was for only a few 

minutes, and is not a problem to the crew. On the positive side 
this controHer was able to quickly reduce the magnitude of the ' 
oscillations and rapidly achieved a steady mass flow rate. 
Comparing these results to our previously listed criteria this 

control ler StantS W3S elected as best for use in the classic PID 
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CLASSIC RESPONSE 

CREW EXCHANGE (25, I, 100) 



Figure 4.2-5 system Response with Weighting (.04,1, .oi) 


InLres^ n9 I t n waS e dl^rerto e con?i™ n ?h t ? i ?h eXer0ise Kas als ° ° 

g a problem that^the ^system was not 


experiencing a problem with that the system was 

sorbent beds require that the e ®J-^ al J- oadlng in the beds. The 
that they don • t^mply continue tnT? °?? le rem ° Ves enou ^ C 
saturation. Figure 4 2 -6shoSo%^ i°; d , unt ^ they reach a 
percentage loading of both^eds olottfd 1 ^ 1119 wi th the 

graph it is easy ?o see that the^fn Versus time * From the 
residual loading problems!^ h b d are not suffering any 


and maintaining desirable cabin conditi 10 regulatin 9 the system 
constants on the response of IL The , effect of the 

air of credibility to the model as fleeted lending an 
be understood that the controllers tested^ 1 ^ - Again let 
of a capable and satisfactorv contJ^f* d chosen in search 

formal optimization study Y controller ' not the result of a 




Figure 4 


2-6 Bed Loading Curves 


4.3 Expert Control Results 

description 


the exact 

Classical Control Results t-k ln ^ ec ^- lon 4.2 called 

modify the expert controller to ^ovidt****** again att empt to 
to lessen the wear on the fan and S some degree of dampening 
the sorbent beds. The understood J-2? < ? r ^ vin ^ the air through 
dampening is that the system mutt £ e ® t f a ^ nt on maximizing 9 
at approximately the 0.0667 kPa set point? ^ C ° 2 levels 

percentage a ^ief, a 

multiplied by a weighting factor to^n * ° han ^ e * Th is belief is 

a new mass flow 

Will generate a high frequency 
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m=m ± ( Jc 1 ( 2 * ^Belief- 1) ) . 


(7) 


controller. The inverse of this is that a small weighting factor 
will result in a lower frequency controller. 

The original controller was designed with kl equal to 0.05. The 
result of this controller when tested with the crew exchange 
scenario is given as Figure 4.3-1. The upper curve corresponds 
to the right hand axis and displays the partial pressure of C0 2 
in the cabin in kPa. The left hand axis goes with the lower 
curve to show the mass flow rate in kg/sec. 


EXPERT RESPONSE 

CREW EXCHANGE (0.05) 



Figure 4.3-1 Dynamic Response with Weighting (0.05) 


The controller exhibits no apparent dampening, and so does not 
appear very suitable for our application. The next course of 
action was to remember, as with our work on the PID controller 
that a lower frequency controller provided smoother mass flow 
rates and an increase in dampening. Following that hunch, the 
value of Kl was lowered to 0.02 and the test was run again. 

The results for the test at kl = 0.02 are given in Figure 4.3-2. 
There is still no evidence of dampening, and the only major 
deviation between the two runs was the fact that the second 
controller was not able to keep the partial pressure of C0 2 as 
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b« & 0 a k ?e r d C u°^a r ?^?a r ct 0 ^ ere£0re 

EXPERT RESPONSE 

CREW EXCHANGE (0.02) 



6 a 

THE (HOURS) 


0.04 


Figure 4.3-2 Dynamic Response with Weighting (0.02) 


The nejct trial was conducted with an even smaller value assigned 

wav to 0 005 th Th? SSt thS w f lghting factor was reduced all the 
y to 0 • 005 . This served to slow the controller's time of 

result^fo^thi^r t0 achi ® ve a . sl ight dampening effect. The 
results for this run are shown m Figure 4.3-3. The Quickest 

dampening however was limited to the region when CO, vSL the* 

lghest . Th i s trend was similairly observed in the PID 

1017 " h ; n .? he f T e 9“«ncy was Slightly to£ high. This 

t S^ tha T thS " elghting factor is close to the desired 
value, and only needs fine tuning. 
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EXPERT RESPONSE 

CREW EXCHANGE (0.005) 



Figure 4.3-3 Dynamic Response with Weighting (0.005) 


The next K1 was then given a value of 0.001 and this test data 
was subjected to the same situation. The graph portraying these 
test results can in Figure 4.3-4. Here the dampening that was 
desired becomes apparent at every level of C0 2 generation. This 
controller however has one major handicap. It was too sluggish 
to appropriately react to the transient cases. The peak at four 
hours and the dip at eight hours both represent undesirable 
deviations from the 0.0667 kPa set point. These deviations are 
short lived, and do not represent a problem for the human 
occupants. The net result being that this is an acceptable 
controller as it meets the basic criteria. 

The final variation on the expert system weighting factor was to 
set k2 - to 0.002. The graph in Figure 4.3-5 represents the 
results of that test. It can be seen that the increase in 
controller frequency enabled the controller to decrease the 
amplitude of the transient spikes. That reduction coupled with 
the fact that the dampening was even more successful made the 
weighting factor of 0.002 appear to be the most capable option 
for the expert controller. The sorbent bed loading curves for 
the test run at kl = 0.002 are also included in Figure 4.3-6. 
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Figure 4 


3-4 Dynamic Response with Weighting 


(O.ooi) 



Figure 4.3-5 Dynamic 


Response with Weighting 


( 0 . 002 ) 
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PPC02 CABIN (KPA) 




EXPERT RESPONSE 



Figure 4.3-6 Bed Loading 


the Classic control ler t i s S not S S D tim^ *5® eX ? ert controller like 

best choice from among several oDtionl* 3 * A ?; t ^ ou 9 h the apparent 
controller for the system. merely a functional and capable 


4.4 Dynamic Case Studies 


DESCRIPTION 

controllers f ove . ex ample S/ the 

Siasarsr S iSB£hkFi s ^ ther 

a heavier load on the same sorlent l^l TMs woul d T .° raated 

frequency? t °The e r esults S for S this e test l3 can 9 be X found 

pressure^on " the * right^axis ^nTthTLsl^fl ^^P^ial 

curve scaled along^he left hand axis l0W ® 13 the bottom 
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PPC02 CABIN (KPA) 



The next case was conducted to determine the natural frequency of 
the controllers. By imparting a single impulse, in this case a 
short term high C0 2 production spike, it is possible to observe 
the systems natural frequency. The results of this test can be 
seen in Figures 4.4-3 and 4.4-4. It can be noticed that the 
expert controller has the higher frequency of the two. That does 
not necessarrily imply that the expert controller has the faster 
response capability, only that it cycles as a higher rate. Also 
in this scenario it is easy to observe the dampening abilities of 
the control systems as they reduce the oscillations amplitudes. 
The final point of interest is the visibility of the half hour 
frequency imparted due to bed switching. It is what is 
responsible for the steady state oscillations visible in the 
graphs . 


CLASSIC RESPONSE 

IMPULSE (25. 1. 100) 



Figure 4.4-3 Classical Response to an Impulse 
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EXPERT RESPONSE 

IMPULSE (0.002) 



Figure 4.4-4 Expert Response to an impulse 

The final scenario examined was the controllers' abilty to handle 
a massive C0 2 production rate. This would simulate a fire in a 
X a £ e S m? tl0n module ' or possibly a leak in the C0 2 accumulator 
. a *? k * T Sh re f ult f of this trial are given in Figures 4.4-5 and 
the'Jwo hS cl ^! slcal astern was able to respond the quicker of 
he two as evidenced by its more rapid increase of the mass flow 

co ^»T+£fi 1 response of the expert system results in the 

CO, partial pressure reaching a value of 14 kPa as opposed to the 
PIO's peak value of 12 kPa. The major consideration however is 
how long before the C0 2 level returns to acceptable limits and 

hoirs^atSr?^ 011613 b ° th Sh ° W th * situation under control by 2 
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CLASSIC RESPONSE 

FIRE LOADING EXCHANGE (25, 1, 100) 



0 2 4 6 8 10 12 


THE (HOURS) 


Figure 4.4-5 Classical Response to a Fire 


EXPERT RESPONSE 

FIRE LOADING (0.002) 



0 2 4 6 8 10 12 


THE (HOURS) 


Figure 4.4-6 Expert Response to a Fire 
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4.5 Conclusions and Recommendations 

CONCLUSIONS 

The first conclusion that can be gathered from this report is 
that the simulation presented is a success. The physical 
phenomena modeled are accurate and respond correctly to parameter 
changes. This implies that the simulation is capable of being 
used as a test bed for evaluating almost any parameter's 
influence on the systems behavior. It is possible to determine 
the effects of the possible disasters (such as a fire), or to 
merely examine how the system operates under normal conditions. 

Both controllers were found to be capable of handling the tasks 
assigned. There is currently no way to evaluate the controllers 
as far as superior capability. Neither was formally optimized, 
and so the limit of their abilities is still not known. 

RECOMMENDATIONS 

It is recommended that a formal optimization of the controllers 
be done. Once optimization is completed, a rigid and weighted 
set of criteria should be drafted. After testing the controls 
with the simulation code, the control schemes could be scored 
against the criteria. Once this is completed, the better control 
system should be implemented as the control scheme of choice. 

Note that a single type of control may not necessarily be the 
best choice. Rather, a control heirarchy where an expert system 
oversees a series of classical controls (or vice versa) might be 
the most effective choice. 
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5.0 APPENDICES 


5.1 Modeling Lecture Summary 


Dr. Byron Jones, professor of mechanical engineering at Kansas 
State University, introduced the design team students to the 
concepts of mathematical modeling and controls. The following is 
an outline of the lecture. 


A. Definitions 

1- Parameters . Parameters which are inherent to the system or 
which are determined from outside the system. Parameters do not 
change state variables but will affect how rate variables relate to 
state variables. 

2. Inputs . Inputs are variables which change state variables. 
The values of Rate Variables are determined from outside the 
system . 

3 . System Relationships. The relationships that describe how the 
values of the rate variables are determined from the state 
variables. (This is the tough part.) 

4. State Variables. Those variables required to define the state 
of a system. Specifying all of the state variables completely 
defines the state of a dynamic system. State variables do not 
change instantaneously. 

5. Rate Variables. Rate variables are those which change state 
variables. The values of rate variables are determined by state 
variables. 


B. Dynamic Systems Simulation 

Development of a set of equations (or other model) that 
describes how a system behaves over time in response to various 
inputs to the system. 


C. Dynamics Systems Analysis 

Use of various techniques to study the dynamic behavior of a 
system (stability, speed of response, etc...) 


D. System of Equations 
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A set of equations must be developed or applied to describe 
the behavior of the system being modeled, and they have the form 


dX 1 

dt 


=F 1 (X 1 I X 2 , 


■ ' X n > y. 




(i) 


E. Solving the Equations 


1 . 

Initial Conditions. 

2. 

Calculate 

Derivatives. 

3 . 

Intearate 

One Time Step. 

4. 

Calculate 

Derivatives. 

5. 

Intearate 

One Time Step. 


F. Example of a Water Supply System 

1. Verbal description of Components 


a) 

State variables 



L- water level (ft) 



N- valve position (number 

of turns open) 

b) 

Inputs 



F 0 - rate at which water is 

demanded (cfm) 

c) 

Parameters 



Lj- float level (10 ft) 

W m - maximum valve speed (10 rpm) 

N m - turns required to fully open valve (20 revolutions) 
F m - maximum water supply flow (100 cfm) 

A - tank area (100 ft 2 ) 

2 . System equations 


dL = £i_£o 
dt A A 


( 2 ) 



(3) 


F^F^N. 

if L < L s and N < N m then 


( 4 ) 


60 



w=w m . 


if L > L. and N > 0 then 


w=-w m , 


otherwise 


W= 0 . 


3 . Initial Approximation 


F i = — xF m 
1 m 


(5) 

( 6 ) 
(7) 


( 8 ) 


4 . Get rid of switch and replace with a P-D controller 

a) Proportional controller 

Output=K 1 x(L set -L) (9) 

b) Derivative controller 

Output=K 2 x-^ (10) 


c) Total output 

Output total =K iX (L seC ~L) + K 2 x^_ 


( 11 ) 
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5.2 Expert Systems Lecture Summary 


Harold Kraus, expert systems consultant at Kansas State University, 
gave a two hour lecture on expert systems. Harold developed an 
expert control model for the 1990-1991 NASA/USRA team. This model 
was developed using CLIPS, which has highly influenced his design 
and the design that is being implemented this year. The one 
problem with using expert systems however in this project is that 
it becomes increasingly difficult to model an Environmental Control 
System when the simulation is written in FORTRAN and the expert 
control is written in CLIPS. CLIPS is a closed system in that it 
cannot communicate with other programming languages except through 
files. Therefore, there comes a problem with passing variables 
between simulation and control modules. Not only is this difficult 
and time consuming to do, it causes the overall model to run slower 
because of the I/O between the disk and memory. To avoid this, 
Harold Kraus implemented the entire model within the CLIPS 
environment. This also made it difficult because anyone who has 
tried to implement a linear project in an expert systems language 
quickly learns that rules do not fire in the order that is 
expected . 

In his lecture, Harold talks about the definitions, properties, 
applications, uses, and examples of an expert system. The 
following is an outline of the lecture: 


I. Expert Systems 

A. Definition of an expert system: 

An expert system is a system that contains expertise 
about knowledge. It makes decisions based on giv en 
information. The separations of the expert system lie in 
the knowledge base (facts and rules) and the engine (CLIPS) 
itself. 

B. Properties of an expert system: 

Properties of an expert system should be consistency and 
reliability. 

C. Applications and uses of an expert system: 

Applications of expert systems include classifications, 
diagnosis, design, and control. Expert systems should be 
used when there are many inputs and outputs and the 
relationship between those inputs and output are inexact 
or incomplete . 

D. Examples of an expert system: 

One of the examples was the discussion of a fuzzy 
controller. Using a fuzzy controller gave a good example 
of a reasoning model and inexact reasoning. The model 
that was used was a simple float valve in a tank of 
water. The following is his example. 
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The initial conditions, as well as time step, run time limit, and 
model parameters are provided to the expert system through an 
initialization data file. Let the file define the following 
values : 


TIME_STEP 

WATER_LEVEL (position) 
VALVE_POS ITION 

At t = 0 minutes, 

LEVEL_ERROR 
WATER_LEVEL (rate) 

Applying the membership functions 
non-zero confidence factor: 


= 0.2 minutes 
= 10 feet 
= 0 revolutions 


= 0 feet 
= 0 ft/min. 

yields just one statement with a 


LEVEL_ERROR is SMALL confidence = 1. 

Applying the operational rules yields the conclusion 

MAINTAIN_VALVE_POS ITION confidence = 1. 

The centroid of MAINTAIN_VALVE_POSITION is zero so 

VALVE_S PEED = 0 rpm. 

Stepping the model through one time step using the calculated valve 
speed changes the input variables such as 

at t = .2 minutes, 

LEVEL_ERROR =-01 feet 

WATER_LEVEL = - 0 .‘ 5 ft/min. 

Applying the membership functions yields the following statements: 

LEVEL_ERROR is POSITIVE confidence = 0, 

LEVEL_ERROR is SMALL confidence = 0.8, 

LEVEL_ERROR is NEGATIVE confidence = 0.2, 

WATER_LEVEL is FALLING confidence = 1, 

WATER_LEVEL is RISING confidence = 0. 

Notice how the confidence values add to 1. They are normalized so 
that we may obtain a percentage weight that gives a numeric 
confidence or belief that a fact is true. The higher the 
confidence number, the more we believe that the fact is true. The 
lower the confidence number, the less we believe that the fact is 
true. 


Applying the operational rules yields the following conclusions: 
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confidence = l, 
confidence = 0, 
confidence = 0.8. 


OPEN_VALVE 
CLOSE_VALVE 

MAINTAIN_VALVE_POSITION 

What could be done in this case is to compare the confidence values 
for OPEN_VALVE , MAINTAIN_VALVE_POSITION, and CLOSE_VALVE. The one 
with a higher confidence value would take precedence. 


II. 1990-91 expert systems project 

A. Explain what last years expert system controlled. 

Last years problem dealt with the Oxygen Generation 
Assembly (OGA) . 

1. The sensory inputs into the Expert System included: 

a) C 2 0 

b) H 2 0 
® ) 0 2 • 

2. The command outputs included: 

a) Reduction of the amount of carbon dioxide 

b) Removal of carbon dioxide from the cabin. 

B. Demonstrate the techniques used in a simple model 

1. Triangular functions 

A triangular function if one method of 
determining percent belief in a fact. With a 
triangular function, one determines the three 
points of the triangle by getting knowledge 
from an expert on the subject matter. 

Triangular functions are useful when trying to 
pull values within a predefined range. An 
example would be to keep the pressure inside a 
tank at a pressure P plus or minus some delta 
P. 

2. Trapezoidal functions 

A trapezoidal function is the same thing as a 
triangular function with the exception that it 
contains four points rather than three. This 
causes the function to take the form of a 
trapezoid with the top sometimes referred to 
as the plateau. 
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5.3 Homework 1 Summary 


INTRODUCTION 

The following is a summary of the first homework assignment 
presented to the advanced design team by Dr. Byron Jones, professor 
of Mechanical Engineering at Kansas State University. It's purpose 
was to introduce the design team students to the concepts of 
mathematical modeling. 


PROJECT DESCRIPTION 

Carbon dioxide is generated by a process at low pressure. A vane 
pump compresses the C0 2 into a storage tank where it is withdrawn 
periodically for use elsewhere. The pump runs at a constant speed 
(1000 rpm) . 3 The vane pump is volumetric, or it pumps a fixed 
volume (40 cm ) of C0 2 into the tank each revolution. The tank is 
perfectly insulated so there is no heat loss to the ambient 
environment. Maximum allowable tank pressure is 300 kPa. 


INITIAL CONDITIONS 

Case 1 . There is no usage of the stored CO z and the inlet 
temperature and pressure remain constant at 40 °C and 2 5 kPa 
respectively. 

Case 2 . Same as case 1 except the inlet pressure decreases with 
time according to the relation below 

P=P xe ( ’ t/a) . (12) 


Stop the simulation after 60 minutes if the limits above are not 
reached prior to that time. 

The symbols used in the above equation are defined as follows. 

P = pressure as a function of time, 

P Q = initial pressure (25 kPa) , 
t = time from start up (minutes) , 
a = a constant (10 minutes). 

Case 3 . Same as case 2 except now C0 2 is used from the reservoir 
at the rate of 0.3 g/s starting at a time 2 minutes and continuing 
for 5 minutes. 
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ASSIGNMENT 

pump system^o/alT Viree'^cases? m ° del ° f the behav ior of the tank 


DYNAMIC MODEL 

1 * Ve rbal Description o f comnnnpn^ 

a) State Variables 

P - Tank Pressure (kPa) 

T - Tank Temperature (^C) 

b) Inputs 

M out ~ rate that C 0 2 is demanded 

Pi - inlet pressure to the pump as a function of time 

c) Parameters 


W„ 


- speed of pump (1000 rpm) 

- displacement of volumetric pump (40 cm 3 ) 


* . 77 — vuiumetric pump (40 

o, tank “ initial pressure of tank (25 kPa) 

^.tank - initial temperature of tank(25°C) 
v t “ Volume of tank (.25 m 3 ) 

Pi, initial - initial pressure at inlet of pump (25 kPa) 
i, initial ~ initial temperature at inlet of pump (40°C) 


System Equations 


*4ank(£) ( t ) X (T^ ( t) -273) 


( 1 ) 


( S *d£) -M^ (£) . xdt 

at 


( 2 ) 


at 


( 3 ) 


3; a nk(t+dt)=- 


^ank(t-dt) 


^tank ( b +dt) XC 


+273 


( 4 ) 
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( 5 ) 


P tank (t + dt) = 


RxT tank (t + dt)xM tank (t + dt) 


^tank 


3 . Simulation (Computer Program) 

PROGRAM PUMPTANK 

REAL MPIN , MPUMP , MTANK , MUSED 

C ... INITIALIZE VARIABLES 
SPEED= 1000. 

RCO2=0 . 1889 
VPUMP=0. 00004 
VTANK=0 .25 
CVC02=0. 657 
CPC02 =CVC02 +RC02 
KC02=CPC02/ CVC02 
H0=RCO2*273 
TIME=0 . 

TTANK=2 5+273 
PTANK=2 5 . 

PPIN=25 . 

TPIN=40 . +273 . 

MTANK= VTANK* PTANK/ (RC02*TTANK) 

UTANK= (TTANK-273 ) *CVC02*MTANK 
MUSED=0 . 0 
DTIME= . 1 

TTANK=UTANK/ (MTANK*CVC02 ) +273 
PTANK=RCO 2 * TT ANK * MTANK/ VTANK 
OPEN (UNIT=11 , FILE= ' HW1 . OUT ' ) 

WRITE (11,*) 'TIME (MIN) TEMP (C) PRESS. (kPa) MASS, TANK' 

C . . . CALCULATE CURRENT STATES 
TO=TIME 

1 IF ( (TIME-TO+ .05) . GT. 1 . ) THEN 

WRI TE ( 1 1 , 1 0 0 ) TIME , TTANK- 2 7 3 , PTANK , MTANK 
TO=TIME 
END IF 

C . . . UNCOMMENT FOR CASE 3 

I F (TIME . GT . 1 . 9 5 . AND . TIME . LT . 7 . 0 5 ) THEN 
MUSED= .0003 
ELSE 

MUSED= 0. 

END IF 

WRITE (*,*) TIME, MUSED 

C ... 

MPUMP=SPEED*VPUMP*PPIN/ (RC02*TPIN) 

TPUMP=TPIN * (PTANK/PPIN) ** (1-1/KC02) 

HIN=MPUMP* ( (TPUMP-273) *CPCO2+H0) 
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HOUT=MUSED* ( (TTANK-273 ) *CPC02+H0) 
TIME=TIME+DTIME 

MTANK=MTANK+ (MPUMP-MUSED) *DTIME 
UTANK=UTANK+ (HIN-HOUT) *DTIME 
TTANK=UTANK/ (MTANK*CVC02 ) +2 7 3 
PTANK=RC02 *TTANK*MTANK/VTANK 

C . . . UNCOMMENT FOR CASE 2 AND CASE 3 
PPIN=2 5 . *EXP ( -TIME/ 1 0 ) 

IF (TIME .GT. 60) GOTO 2 
C 

IF(PTANK. LT. 300. .AND.TTANK.LT. (250+273) )GOTO 1 
2 WRITE (11, 100) TIME, TTANK-273, PTANK,MTANK 

100 FORMAT (F5 . 2 , 3 ( 2X , F10 . 4 ) ) 

END 


4 . Results 

Case 1 . The conditions for the first case described a constant 
inlet temperature and pressure. Figure 5.3-1 shows results of the 
first case. Notice that because a constant inlet pressure is 
maintained, the pressure within the tank increases linearly. As 
the pressure increased with time, the temperature seemed to 
approach an asymptotical value of about 394 K. The simulation was 
halted after 53 minutes due to the pressure exceeding its limit. 



Temperature Pressure 


Figure 5.3-1: Homework 1 Case 1 Data 
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Case 2. Case 2 was the same as case 1 except for a decreasing 
inlet pressure according to a given function. As shown in Figure 
5.3-2, the temperature and pressure seemed to level out at 360 K 
and 77 kPa respectively. The pressure within the tank responded as 
would be anticipated with an exponentially decreasing input 
function for the inlet pressure. The temperature increased much 
more rapidly for case 2 from 0 to 10 minutes, and it. leveled out. at 
a values about 30 K less than the previous simulation. No limits 
on the tank were exceeded so the simulation automatically halted at 
60 minutes. 



Temperature Pressure 


Figure 5.3-2: Homework 1 Case 2 Data 


Case 3 . Case 3 added a withdrawal of .3 g/s from the tank from 2 
minutes to 7 minutes. The results from this simulation, shown in 
Figure 5.3-3, were similar to that of Case 2, and the effects of 
the new condition seemed negligible. If figures for Cases 2 and 3 
were examined closely from 2 minutes to 7 minutes,, slightly less 
slope would be seen in both responses due to the drain on the tank. 
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Temperature (K) 



Temperature Pressure 


Figure 5.3-3: Homework l Case 3 Data 
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Pressure (kPa) 



5.4 Homework 2 Summary 


INTRODUCTION 

The following is a summary of the second homework assignment 
presented to the advanced design team. It's purpose was to 
introduce the design team students to the concepts of simulation 
diagrams and numerical integration of differential state equations 
using methods other than the Euler time step. 


PROJECT DESCRIPTION 

The system examined is shown in Figure 5.4-1. It is a simplified 
automobile suspension consisting of h of the automobile mass, M x , 
the suspension system mass, M 2 , the suspension spring, K 1( the 
elasticity of the tire, K 2 , and the shock absorber, B. The 
vertical displacement of the road is input to the system by a 
force, f(t), acting on the suspension mass. A major simplification 
is made by assuming the tire to never leave the road surface. The 
coefficients and forcing function are given as follows. 

= 250 kg K x = .55 kg/s 2 

M 2 = 30 kg K 2 = 20 kg/s 2 

B = 1.5 kg/s f (t) = 10 e 1 sin (27rt/300) N 



x (t) 
l 


x (t) 
2 


Figure 5.4-1: Simplified Automobile Suspension System 
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ASSIGNMENT 


1. Write dynamic equations that simulate this system. Summing 
forces on each mass should result in two second-order 
differential equations. 

2. Draw a simulation diagram from these simulation equations. 
How many integrators should it have? 

3. Write a set of state equations describing this system. 

4. Produce a set of four plots: M x position, M z position, Mj 

velocity, and M 2 velocity. To produce these plots, write or 
use a computer program to perform the simulation. 


SOLUTION 

1. Summing forces on M x and M 2 gives equations 1 and 2 


d 2 Xi 

M, r = -0 

dt 2 


dx, dx„ 


dt dt 


1(— l 2 ) 


d 2 x 

M 2 f = f (t) - B| 

dt 2 


dx. dx. 


dt dt 


i (— 2 i ) ~— 2—2 


( 1 ) . 


( 2 ) 


2. Figure 4.5-2 is the simulation diagram with four integrators. 



Figure 5.4-2: Homework 2 Simulation Diagram 
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3 . 


The following 3-6 are a set of state equations obtained 
directly from the simulation diagram. 


= ^(-Bx 1 + Bx 2 + x 3 ), 


Mi 


X 2 - -^(b^-Bx, xj. 


m 2 


*3 = -SA+IA. 


* 4 = f (t) -K^ + K^-K^ 


( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 


Another set of state equations is arrived at by first making 
the definitions 




= 


—2 


= X, 


( 7 ) 


( 8 ) 


Substituting these into the second-order differential 
Equations 1 and 2 gives 


—3 



( 9 ) 


—4 





— 2—2 


( 10 ) 


4. Plots of M : position, M 2 position, M x velocity, and M 2 velocity 
can be seen in Figures 5.4-3 through 5.4-6. The solution was 
arrived at by solving the equations every .1 seconds for 300 
seconds using the 4 <h /5 th order Runga-Kutta-Fehlburg method 
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(RKF45 subroutine) . The results are underdamped solutions. 
The "bump” causes the suspension M 2 to vibrate for roughly 4 
minutes, while the automobile Mi has only lost half of its 
vibrational amplitude by the end of the 5 minute simulation. 



Figure 5.4-3: M x Position 



Figure 5.4-4: M 2 Position 
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Velocity (m/s) 


200 


250 


3 


iO 

Time (s) 

Figure 5.4-5: ^ Velocity 



Figure 5.4-6: M 2 Velocity 


300 


5 





noon 


Computer Program 


C$INCLUDE DEQ 
C$INCLUDE RKF45 

PROGRAM HW1 
INTEGER IWORK ( 5 ) 

DOUBLE PRECISION X ( 4 ) , T , TOUT , RELERR, ABSERR, WORK ( 27 ) 
EXTERNAL DEQ 


N=4 

T=TOUT=0 . 

X ( 1 ) =X ( 2 ) =X ( 3 ) =X ( 4 ) =0 . 

IFLAG=1 
RELERR=lE-08 
ABSERR=0 . 

OPEN (UNIT=2 , FILE= ' HW1 . DAT ' ) 

100 FORMAT (F6.1,4E12.4,I3) 

WRITE (2, 100) T,X(1) ,X(2) , X ( 3 ) ,X(4) , I FLAG 
D010I=1, 3000 
TOUT=I/10 . 

5 CALL RKF4 5 ( DEQ, N,X,T, TOUT, RELERR, ABSERR, I FLAG, WORK, IWORK) 

IF (IFLAG. NE . 2 ) THEN 

WRITE ( 6 , ' (A, 12 , A, F6 . 1 , A) ') ' IFLAG= IFLAG, ' at',T, 

& ' seconds. . .trying again' 

IF (IFLAG. EQ. 7) IFLAG=2 
GOTO 5 
ENDIF 

WRITE ( 2 , 100) T , X ( 1 ) ,X(2) ,X(1) ,X(2) , IFLAG 
10 CONTINUE 

END 

SUBROUTINE DEQ (T , X, XDOT) 

DOUBLE PRECISION T, X ( *) , XDOT ( *) , F 
REAL B , K1 , K2 , Ml , M2 


B=1 . 5 
Kl= . 55 
K2=2 0 
Ml=250 
M2=3 0 

F=10*EXP(-T) *SIN(T*2*ACOS (-1. )/300) 

XDOT ( 1 ) =X ( 3 ) 

XDOT ( 2 ) =X ( 4 ) 

XDOT ( 3 ) = (B* (X (4 ) -X ( 3 ) ) +K1* (X (2 ) -X ( 1) ) )/Ml 
XDOT ( 4 ) = ( B* ( X ( 3 ) -X ( 4 ) ) +K1* ( X ( 1 ) -X ( 2 ) ) -K2 *X ( 2 ) +F ) /M2 
XDOT ( 1) = (B* (X (2 ) -X ( 1) ) +X (3 ) ) /Ml 
XDOT ( 2 ) = ( B* (X(l)-X(2)) +X ( 4 ) ) /M2 
XDOT ( 3 ) =K1* (X (2 ) -X ( 1) ) 

XDOT ( 4 ) =K1* (X(l) -X(2) ) -K2*X(2)+F 
RETURN 
END 
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5.5 Classical Controls Source Code 


C$NOEXT 

C$NOWARN 

C*************************************************************** 


C PROGRAM NAME : 

C GROUP NAME : 

C CREATED BY : 

C REVISED: 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C CALLS : 

C CALLED FROM: 

C OUTPUT FILE: 

C 


c 

/I 

LIST OF 

VARIABLES : 



c 

c 

INT 

TIME = CURRENT TIME 

c 

REAL 

MDOT = MASS 

FLOW OF AIR INTO SYSTEM IN Kg/SEC 

c 

REAL 

PC02 = PARTIAL PRESSURE OF C02 IN CABIN AIR 

c 

INT 

BLOSPEED = 

: BLOWER SPEED, ONE OF THREE 

c 


1 

-- 

- NORMAL 

c 


2 

-- 

- DEGRADED 

c 


3 

-• 

- EMERGENCY 

c 

INT 

BEDSWITCH 

= 

TIME INTERVAL FOR SWITCHING SORBENT BEDS 

c 




— SEC 

c 

REAL 

PTANK = C02 

ACCUMULATOR TANK PRESSURE — KPa 

c 

REAL 

PUMPS PE ED 

= 

C02 ACCUMULATOR TANK PUMP SPEED — RPM 

c 

REAL 

PTANK SET 

= 

DESIRED PRESSURE OF C02 TANK — KPa 

c 

REAL 

PTANK_OLD 

= 

C02 TANK PRESSURE FROM LAST TIME STEP 

c 




— KPa 

c 

INT 

SYS 

= 

SYSTEM STATUS, 1 = ON, 0 = OFF 

c 

INT 

SYSTIME 

= 

TIME SINCE SYSTEM WAS LAST TURNED ON, 

c 




VARIES FROM 0-1800 

c 

REAL 

CO 2 REM 

= 

CO 2 REMOVED BY SYSTEM 

c 




IF VALVE = 0, C02REM = MC02R 

c 




IF VALVE = 1, C02REM = 0 

c 

REAL 

PB 

= 

PRESS OF DESORBING SORBENT BED — KPa 

c 

REAL 

PBC02 

= 

EQUILIBRIUM PRESS OF C02 IN DESORBING 

c 




SORBENT BED — KPa 

c 

REAL 

MAIR 

= 

MASS OF AIR LEAVING DESORBING SORBENT 

c 




BED — Kg/SEC 

c 

REAL 

MC02R 

= 

MASS OF C02 LEAVING DESORBING SORBENT 


C BED — Kg/SEC 


CCONTROL. FOR 
CLASSICAL CONTROLS 
STAN KROEKER 
02MAR92 - CREATED 

04MAR92 - COMPLETED FIRST TIME 

16MAR92 - ADDED SYSTEM CONTROL AND CABIN 

SUB. 

17MAR92 - ADDED VALVE BETWEEN PUMP AND 

TANK 

22MAR92 - ADDED BOSCH REACTOR 

24MAR92 - TESTED ALL CONTROLS AND REMOVED 

CABIN. FOR - SATISFIED WITH PROGRAM 


READS CDATA 
C. FOR 

WRITES TO CDATA 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


REAL MS BOUT = MASS OF AIR LEAVING ADSORBING SORBENT 

BED — Kg/SEC 

INT VALVE = STATE OF VALVE BETWEEN PUMP AND TANK 

1 DIRECTS GAS TO CABIN 
0 DIRECTS GAS TO TANK 

REAL MTANKUSED = MASS OF C02 USED BY BOSCH REACTOR Kg/SEC 
DESCRIPTION: 

PROGRAM READS FROM CDATA TO DETERMINT THE CURRENT STATUS OF 
THE CONTROL PARAMETERS, CALCULATIONS ARE THEN MADE TO CHANGE 
THE CONTROL VARIABLE IF NEEDED. 

TO CONTROL THE PARTIAL PRESSURE OF C02 IN THE CABIN: 

TURNS THE SYSTEM ON FOR 1800 SEC IF THE PC02 > .4KPa 
TO CONTROL THE C02 ACCUMULATOR TANK PRESSURE: 

CHANGE THE C02 ACCUMULATOR PUMP SPEED 


C* ************************************************************** 


SUBROUTINE PID (SIMTIME , DT , MDOT , PRESS , TEMP , PHI , TIME , PC02 , 
& BLOSPEED, SYSTIME2 , PTANK, PUMPS PEED, PTANK_SET, 

& PTANK_OLD , S YS , SYSTIME , CO 2 REM , PB , PBC02 , MC02R, 

& MC02P, VALVE, MTANKUSED) 

REAL PTANK, PUMPS PEED, PTANK_SET , PTANK_OLD, PRESS 
DOUBLE PRECISION MDOT, MC02R, MC02P 
REAL PBC02 

DOUBLE PRECISION MTANKUSED, C02REM 

INTEGER TIME , BLOSPEED , BEDSWITCH , SIMTIME , DT , SYS , SYSTIME 
INTEGER VALVE, SYSTIME2 

C ... PID CONTROL OF C02 REMOVAL SYSTEM, BASED ON THE PARTIAL 
C PRESSURE CO 2 IN THE CABIN. 

SYS = 1 

IF (TIME . EQ . 0 ) PC02_0LD=PC02 
PC02_SET = 0.0667 
K1 = 25. 

K2 = 1. 

K3 = 100. 

ERR = PC02 -PC02_SET 
DERR = PC02 -PC02_0LD 
ERRDER = DERR/DT 

ERRINT=DT* (ERR+ (PC02_0LD-PC02_SET) )/2 
DMDOT= (ERR/K1+ERRDER/K2+ERRINT/K3 ) 

MDOT = MDOT + DMDOT 

IF (MDOT. GE. 0.5) MDOT =0.5 
IF (MDOT. LE.1E-3) MDOT = IE-3 
PC02 OLD = PC02 


SYSTIME2=SYSTIME2+DT 

SYSTIME=SYSTIME+DT 
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C . . . CONTROL OF VALVE BETWEEN C02 PUMP AND C02 TANK 
IF ( ( PB-PBC02 ) .GT. ( . 01*PBC02) ) THEN 
PUMPS PEED = 3000 
VALVE = 1 

C02REM = MC02R-MC02P 
GOTO 200 
ELSE 

VALVE = 0 
C02REM = MC02R 
PUMPS PEED = 0 
ENDIF 

200 CONTINUE 

C MASS FLOW DEMAND OF THE BOSCH REACTOR 

IF (PTANK. GE .101.125) THEN 
MTANKUSED = 4.633333E-5 
ELSE 

MTANKUSED = 0. 

ENDIF 

RETURN 

END 
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5.6 Expert Systems Source Code 


EXPERT CONTROLS FOR NASA DESIGN PROJECT 


Created : February 21, 1992 

Programmers : Michael W. Honas 

: Robert A. Swenson 

With help by : Dr. David A. Gustafson 
: Stan Kroeker 

Code : '/nasa/mike/expert . clp 

Source Machine : DEPOT 

Modified : April 3, 1992 -removed unused variables and 

modified the controls 

combined VALVE-TO-TANK & VALVE-TO-CABIN 
into VALVE-TO-CABIN 


Description: 

This will monitor the simulation and write 
new control parameters to the 'rates' file. The only 
expert controlled device this far are the tank pressure 
and the cabin pressure. Other devices being controlled 
by on/off switches are as follows: C02, valve to the cabin, 
and bosch flow rate. 


TRIANGULAR FUNCTIONS — 

The following fuzzy sets manually define the two triangles that 
will determine if the pressure in the tank/cabin is too high 
too low, or anywhere inbetween. The beauty of this method is 
that these triangles can easily be changed to tweak the results. 
Further more, this can be done on the fly while the program is 
running. Simply edit the file and then change the values and 
save them back out. 


(deffacts start 

(state openedata) 
(fuzzy tank-pressure 
(fuzzy tank-pressure 
(fuzzy cabin-pressure 
(fuzzy cabin-pressure 


) 


low -150 300 400) 
high 300 400 650) 
low -100 .058 .075) 
high .058 .075 100) 
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OPEN-edata 

Read in the values from the Fortran and C simulation modules, 
"edata" is a temporary file that is used to pass variables 
between the simulation and the clips control. This is done 
simply because we cannot link CLIPS in with any other compiler. 
This is one of the few draw backs is that it is almost a closed 
system. 


(defrule open-edata 

?s <- (state openedata) 

=> 


(printout t 
(retract ?s) 
(open "edata 
(assert 
(assert 
(assert 
(assert 
(assert 
(assert 
(assert 
(assert 
(assert 
(assert 
(assert 
(assert 
(assert 
(assert 
(assert 


"Entering open expert data" crlf) 


" data " r+") 
(time-step =( 

(mass-in =( 

(partial-co2 =( 
(tank-pressure = 
(pump-speed =(read 
(sys =(read 

(sys-time =(read 
(sys-2time =(read 
(co2-rem =(read 
(pb 

(pbco2 = 

(mco2p =(read 

(mco2r =(read 

(valve =(read 

(mtankused =(read 


read data) ) ) 
read data) ) ) 
read data) ) ) 
(read data))) 
data) ) ) 
data) ) ) 
data) ) ) 
data) ) ) 
data) ) ) 

(read data))) 
(read data) ) ) 
data) ) ) 
data) ) ) 
data) ) ) 
data) ) ) 


(close data) 

(assert (state control-pp) ) 


CONTROL- PC02 

This rule controls the partial pressure of the carbon dioxide 
and also takes care of the system time. 


(defrule control-pco2 

?s <- (state control-pp) 
?f <- (sys ?sy) 

?h <- (sys-time ?st) 

?i <- (sys-2time ?s2t) 

?j <- (mass-in ?mi) 


81 




(partial-co2 ?pc) 
(time-step ?ts) 


(printout t "Entering CONTROL- PC02" crlf) 
(retract ?s) 

(if (= ?sy 1) then 

(retract ?i) 

(assert (sys-2time =(+ ?s2t ?ts) ) ) 


(if (and (and (= ?sy 1) (>= ?st 1800) ) (<= ?pc .4)) then 

(retract ?f) 

(retract ?h) 

(retract ? j ) 

(assert (sys 0) ) 

(assert (sys-time 0)) 

(assert (mass-in 0)) 

else 

(if (and (= ?sy 0) (>= ?pc .4)) then 

(retract ?f) 

(retract ?j) 

(assert (sys 1)) 

(assert (sys-time 0)) 

(assert (mass-in .5)) 

else 

(if (and (and (= ?sy 1) (>= ?st 1800) ) (> ?pc .4)) then 
(retract ?h) 

(assert (sys-time 0)) 

(assert (state valve) ) 

else 

(if (and (= ?sy 1) (< ?st 1800)) then 
(retract ?h) 

(assert (sys-time =(+ ?ts ?st) ) ) 

)))) 

(assert (state valve)) 

) 


VALVE-TO-CABIN 

This system will turn the valve on and off and assert a 
pump speed of 3000 if necessary. 


(defrule valve-to-cabin 
?s <- (state valve) 
(pb ?pb) 
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(pbco2 ?p2 ) 

?f <- (pump-speed ?ps) 
?g <- (valve ?va) 

?h <- (co2-rem ?co) 
(mco2p ?mcp) 

(mco2r ?mcr) 


=> 

; (printout t "Entering VALVE-TO-CABIN" crlf) 
(retract ?s) 

(retract ?g) 

(retract ?h) 

(if (> (- ?pb ?p2 ) (* ?p2 0.01)) ; if pb » p2 

then 

(retract ?f) 

(assert (pump-speed 3000) ) 

(assert (valve 1) ) 

(assert (co2-rem =(- ?mcr ?mcp) ) ) 

(assert (state bosch) ) 

else 

(assert (valve 0)) 

(assert (co2-rem ?mcr) ) 

(assert (state calc)) 

) 

) 


CALC-MEMBERSHIP-FOR-PUMP 

Calculate the percent belief that the tank pressure or cabin 
pressure is low, right or high. 


(defrule calc-membership 
(state calc) 

(fuzzy ?var ?qual ?low ?mid ?high) 
(?var ?val) 


; (printout t "Entering CALC-MEMBERSHIP-FOR-PUMP" crlf) 

(if (and (> ?val ?low) (<= ?val ?mid) ) then 

(assert (member ?var ?qual =(/ (- ?val ?low) (- ?mid ?low) ) ) ) ) 
(if (and (> ?val ?mid) (< ?val ?high) ) then 

(assert (member ?var ?qual =(/ (- ?high ?val) (- ?high ?mid) ) ) ) ) 

) 


RULES 

This rule decides whether an action should take place. 
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That action is based on the belief of better than 50% 


(defrule rules 

(state calc) 

?f <- (member ?variable ?quality ?val&:(> ?val .5)) 

=> 

; (printout t "Entering RULES" crlf) 

(retract ?f) 

(assert (action ?variable ?quality ?val)) 


NO-CHANGES 

Conversly, this rule decides that no action should take place 
if the belief is 50% or less. 


(defrule no-changes 
(state calc) 

?f <- (member ?variable ?quality ?val&: (<= ?val .5)) 

=> 

; (printout t "Entering NO-CHANGES" crlf) 

(retract ?f) 


ACTI ON- S LOW-MOTOR 

This rule will slow the motor by some variable amount 
based on the percent belief. 


(defrule action-slow-motor 
(state calc) 

?fm <- (action tank-pressure high ?val) 

?ps <- (pump-speed ?wps) 

=> 

; (printout t "Entering ACTION-SLOW-MOTOR" crlf) 

(retract ?fm) 

(retract ?ps) 

(if (< ?val (/ ?wps 3000)) then (assert (pump-speed = 

(- ?wps (* 100 ?val)))) 
else (assert (pump-speed 0))) 
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) 


ACT I ON - FAS T-MOTOR 


This rule will speed up the motor by 
based on the percent belief. 


some variable amount 


(defrule action-fast-motor 
(state calc) 

?fm <- (action tank-pressure low ?val) 

?ps <- (pump-speed ?wps) 

; (printout t "Entering ACT I ON - FAS T-MOTOR " crlf) 

(retract ?fm) 1 

(retract ?ps) 

(if (> ?val (/ ?wps 3000)) then (assert (pump-speed = 

. . (+ ?wps (* 100 ?val) ) ) ) 

else (assert (pump-speed 3000))) 


ACTION-S LOW-MDOT 


This rule will slow down the 
based on the percent belief. 


mass flow rate by some variable amount 


(defrule action-slow-mdot 
(state calc) 

?fm <- (action cabin-pressure high ?val) 

?ps <— (mass— in ?mi) 

; (printout t "Entering ACTION-SLOW-MDOT" crlf) 

(retract ?fm) 

(retract ?ps) 

(if (< ?val (/ ?mi .5)) then (assert (mass-in = 

. , (- ?mi (* .05 ?val) ) ) ) 

else (assert (mass-in 0) ) ) 


, . ACTION-FAST-MDOT 

based'o^the'percent “ 0 “ rate by S °" e variabla 
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(defrule action-fast-mdot 
(state calc) 

?fm <- (action cabin-pressure low ?val) 

?ps <- (mass-in ?mi) 

=> 

; (printout t "Entering ACTION-FAST-MDOT" crlf) 

(retract ?fm) 

(retract ?ps) 

(if (> ?val (/ ?mi .5)) then (assert (mass-in = 

(+ ?mi (* .05 ?val) ) ) ) 
else (assert (mass-in .5))) 

) 


CALC-DONE 

This will end the calculations and adjustments made for 
the fuzzy logic sets. Notice that there is a salience of -10. 
This is so that the calculations can complete recursively until 
they are all done. 


(defrule calc-done 

(declare (salience -10) ) 

?s <- (state calc) 

=> 

; (printout t "Entering CALC-DONE" crlf) 
(retract ?s) 

(assert (state bosch) ) 


BOSCH-MASS -FLOW 

This rule will change the mass in the tank used depending 
on whether the tank pressure is above or below 101.125. 


(defrule bosch-mass-flow 

?s <- (state bosch) 
(tank-pressure ?tp) 

?f <- (mtankused ?mt) 


(printout t "Entering BOSCH-MASS-FLOW" crlf) 
(retract ?s) 

(retract ?f) 

(if (>= ?tp 101.125) then 

(assert (mtankused .000046333)) 
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) 

(if (< ?tp 101.125) then 

(assert (mtankused 0) ) 

(assert (state bosch-sys) ) 


BOSCH-SYSTEM 


(defrule bosch-system 

?s <- (state bosch-sys) 
(tank-pressure ?tp) 

?g <- (sys ?sy) 

?h <- (sys-time ?st) 

?i <- (mass-in ?mi) 


(printout t "Entering BOSCH-SYSTEM" crlf) 
(retract ?s) 

(if (and (< ?tp 137.000) (= ?sy 0)) then 

(retract ?g) 

(retract ?h) 

(retract ?i) 

(assert (sys 1)) 

(assert (sys-time 0)) 

(assert (mass-in .5)) 

) — 
(assert (state writeedata) ) 


WRITE-edata 

Write in the values for the Fortran and simulation modules. 


(defrule write-edata 

?s <- (state writeedata) 
?ts <- (time-step ?wts) 

?mi <- (mass-in ?wmi) 

?pc <- (partial-co2 ?wpc) 
?tp <- (tank-pressure ?wtp) 
?ps <- (pump-speed ?wps) 

?sy <- (sys ?wsy) 

?st <- (sys-time ?wst) 
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* st 2 <- (sys-2time ?wst 2 ) 
* co <- (co2-rem ?wco) 

?Pb <- (pb ?wpb) 

?p2 <- (pbco2 ?wp2) 

?xna <- (mco2p ?wxna) 

?mc <- (mco2r ?wmc) 

?va <- (valve ?wva) 

?mt <- (mtankused ?wmt) 


(printout t "Entering write expert data" orlf) 

(retract ?s) 

(retract ?ts) 

(retract ?mi) 

(retract ?pc) 

(retract ?tp) 

(retract ?ps) 

(retract ?sy) 

(retract ?st) 

(retract ?st2) 

(retract ?co) 

(retract ?pb) 

(retract ?p 2 ) 

(retract ?ma) 

(retract ?mc) 

(retract ?va) 

(retract ?mt) 

(system "rm edata") 

(open "edata” data "w") 

(printout data ?wts crlf) 

(printout data ?wmi crlf) 

(printout data ?wpc crlf) 

(printout data ?wtp crlf) 

(printout data ?wps crlf) 

(printout data ?wsy crlf) 

(printout data ?wst crlf) 

(printout data ?wst2 crlf) 

(printout data ?wco crlf) 

(printout data ?wpb crlf) 

(printout data ?wp2 crlf) 

(printout data ?wma crlf) 

(printout data ?wmc crlf) 

(printout data ?wva crlf) 

(printout data ?wmt crlf) 

(close data) 
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5.7 Simulation Source Code 


C$NOEXT 
C$NOWARN 
C$TIME=36000 
C$ INCLUDE SIM 
C$ INCLUDE PID 
C$ INCLUDE EXP 
***** 

* CREATED BY: STAN KROEKER 

* CREATED: APRIL 1, 1992 

* 

* THIS PROGRAM IS TO TEST THE SIMULATION AND THE CONTROLS 

* 

***** 

PROGRAM COOL 

DOUBLE PRECISION MD0T,MAIR,MC02R, MC02P,MTANKUSED, C02REM, MC02 
INTEGER S IMTIME , DT , TIME , BLOSPEED , BEDSWITCH , SYSTIME , S YS , VALVE , I 
INTEGER SYSTIME2 

C . . . REAL DESSICANT BED VALUES 

REAL MDBADS (2,2), MDBTANK ( 2 ) 

C REAL SORBENT BED AND SORBENT PUMP VALUES 

REAL MSBTANK ( 3 ) ,MSBADS(2,3) ,TBED(2,3) 

C . . . REAL C02 TANK VALUES 

DOUBLE PRECISION MTANK , MT IN, MTOUT 

* INITIALIZE VARIABLES 

SIMTIME = 0 
SYSTIME2 = 0 
DT = 6 
MDOT = 0 
PRESS = 101.125 
TEMP =300 
PHI = .78 
TIME = 0 
PC02 = .0667 
BLOSPEED = 0 
BEDSWITCH = 0 
PTANK =300 
PUMPSPEED = 0 
PTANK_SET = 350 
PTANK_OLD = 101.125 
SYS = 0 
SYSTIME = 0 
C02REM = 0 
PB = 101.125 
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o o 


0667 


PBC02 = . 

MAIR = 0 
MC02R = 0 
MC02P = 0 
VALVE = 0 
MTANKUSED = 0 
1=0 
J=0 

CPAIR=1 .006 
MDBTANK ( 1 ) =1 0 0 . 

MDBTANK (2 ) =100 . 

MDBADS (1,1) =0 .001 
MDBADS (1,2) =0.001 
MDBADS (2,1) =0.001 
MDBADS (2, 2) =0.001 

C . . . BLOWER/ PRECOOLER SUBROUTINE VARIABLE INITIALIZATION 

10 CPH20=4 .184 

EPSILON= . 8 
PAIR=1. 1614 
TWATER=288 
VELAIR=3 . 

. . . SORBENT BED AND SORBENT PUMP SUBROUTINE VARIABLE 
. . . INITIALIZATION 

CPAIR= 1.006 
CVAIR=. 719 
CVBED=1 . 

CVC02= . 7 
HC02=572 . 

HEAT=0 . 

MCO2=0 . 

MAIR= . 1 
MS BADS (1,1) =0 . 

MS BADS ( 1 , 2 ) =0 . 

MS BADS ( 1 , 3 ) =0 . 

MS BADS (2,1) =0 . 

MS BADS (2,2) =0 . 

MS BADS (2,3) =0 . 

MSBTANK(l) = 30. 

MSBTANK ( 2 ) = 30. 

MSBTANK ( 3 ) = 30. 

PBED=101 . 325 
PSPIN=1 0 1.325 
RAIR=. 287 
RC02=. 1889 
SPEED=3000 
TBED ( 1 , 1 ) =3 0 0 . 

TBED (1,2)=300. 

TBED(1 , 3 ) =3 00 . 
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TBED(2 , 1)=300. 

TBED (2,2)=300. 

TBED ( 2 , 3 ) =3 00 . 

VBED= . 1 

VPUMP = 0.0006 

C . . . C02 TANK SUBROUTINE VARIABLE INITIALIZATION 


MTANK=2 . 

TTANK=300. 

VTANK=1 . 

MTOUT=0 . 

WRITE (*,*) 'Working. . . ' 

DO 100 TIME = 0, 7500000, DT 

IF (MOD (TIME , 3600) .EQ.0) THEN 

WRITE(*, *) 'The time now is ' , J, ' hours SYS = ' , SYS 

J=J+1 


IF ( J . EQ . 25 ) STOP 


ENDIF 


C 

C DETERMINE CONTROLLER TO UTILIZE 
C 


& 

& 

& 


CALL PID ( S IMTIME , DT , MDOT , PRESS , TEMP , PHI , TIME , PC02 , 

BLOSPEED, SYSTIME2 , PTANK, PUMPSPEED, PTANK_SET , 
PTANK_OLD, SYS , SYSTIME , C02REM, PB, PBC02 , MC02R, 
MC02P , VALVE , MTANKUSED) 


& 

& 

& 


CALL EXP ( S IMTIME , DT , MDOT , PRESS , TEMP , PHI , TIME , PC02 , 

BLOSPEED, SYSTIME2 , PTANK, PUMPSPEED, PTANK_SET , 
PTANK_0LD , SYS , SYSTIME , C02REM , PB , PBC02 , MC02R, 
MC02 P , VALVE , MTANKUSED ) 


C 

C CALL MAIN SIMULATION 
C 

CALL SIM (S IMTIME , DT , MDOT , PRESS , TEMP, PHI , TIME , 

& PC02 , BLOSPEED, MC02 ,MAIR, PTANK, PUMPSPEED, 

& PTANK_SET , PTANK_OLD , SYS , SYSTIME2 , C02REM , 

& PB , PBC02 , MC02R, MC02 P , VALVE , MTANKUSED , CPAIR , 

& MDBTANK , MDBADS , CPH20, EPSILON , PAIR , TWATER , VELAIR , 

& CVAIR , CVBED , CVC02 , HC02 , HEAT , MS BADS , MSBTANK , 

& RAIR , RC02 , TBED , VBED , VPUMP , MTANK , TTANK , VTANK , ADSORB , 

& DESORB, C02 ADSORB, C02 DESORB) 

1 = 1+1 


100 CONTINUE 
END 
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c 

c ... PROGRAM NAME: SIM. FOR 

P NAME: CLASSICAL CONTROL 

c :::.? REATED BY: ™-l m : sn ^er, STA n kr 0E ker 

c ... REVISED: 03/10/92 (PS) 2 

c ... SUBROUTINE CALLS: DESSBED. FOR 

c BLOWCOOL.FOR 

c SORBED. FOR 

c SORPUMP . FOR 

c * " * C02TANK. FOR 

C ... CALLED * FROM i * ‘ CO^REM^FO^ * F ° R ~~ (IN EXPERT SI M ONLY) 

C ... INPUT FILE: CREAD . FOR 

C ... OUTPUT FILE: CWRITE.FOR 

C .*.*.*.*****[ LIST OF VARIABLES 

C ... SEE SUBROUTINES * FOR * LOCAL * VARIABLE * LIST 

^ 

CSNOEXT 

C$NOWARN 

C$INCLUDE DESSBED 
C$ INCLUDE BLOWCOOL 
C$ INCLUDE CWRITE 
C$ INCLUDE SORBED 
C$INCLUDE SORPUMP 
C$ INCLUDE C02TANK 
C$ INCLUDE CABIN 
C$TIME=3600 

^SUBROUTINE ^IMTIME, DT, MDBIN, PDBIN, ™£™IN ; SYSTIME, 
l fMC02R^MC02pf VALVe! MTOOT^PAIR 

& DESORB , C02 ADSORB , C02 DESORB )^ TANK ' TTANK ' VTANK ' ADS0RB / 

c *.. Real DESSICANT BED VALUES 

?? AE MDBADS (2,2), MDBTANK ( 2 ) 

OUBLE PRECISION MH20, MH20R, MDBIN, MDBOUT1 ,MDBOUT2 ,MDBOUT3 
c ••• REAL BLOWER/ precooler values 
DOUBLE PRECISION MBCIN , MBCOUT 
c . . . REAL SORBENT BED AND SORBENT PUMP VALUES 
REAL MSBTANK ( 3 ) , MSBADS (2 3) TBFD f o d 

DOUBLE PRECISION MAIR, MC02 J MC02R,MC02P,MSB0UT,MSP0UT, C02REM 
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C ... REAL C02 TANK VALUES 


DOUBLE PRECISION MTANK , MTIN , MTOUT 

C ... DECLARE ALL INTEGERS FOR THE SIMULATION 

INTEGER DT , T , TIME , BLOSPEED , BEDSWITCH , SIMTIME , SYS , SYSTIME , VALVE 

C . . . RUN CO 2 REMOVAL ASSEMBLY SIMULATION AT A TIME STEP DT FOR A 
C . . . TIME DURATION T 

CALL CABIN ( TDBIN , PDBIN , PHI IN , DT , SYSTIME , C02REM , PC02 ) 

IF (SYS . EQ. 1) THEN 

CALL DESSBED (MDBIN , PDBIN, PHI IN, TDBIN, MDBOUT1 , PDBOUT1 , PHI 
OUT1 , TDBOUT1 , CPAIR, DT , MDBADS , MH 2 O , MH20R, MDB 
TANK , 1 , TIME , ADSORB) 

CALL BLOWCOOL (MDBOUT1 , PDBOUT1 , TDBOUT1 , MBCOUT , PBCOUT , TBCO 
UT , CPH20 , EPSILON , PAIR, TWATER, VELAIR) 

CALL SORBED (MBCOUT , PBCOUT , TBCOUT , MSBOUT , PSBOUT , TSBOUT , CP 
AIR,-CVBED , DT , HC02 , MS BADS , MSBTANK , PC02 
, PVC020UT , TBED, TIME , C02ADS0RB,MC02R) 

CALL DESSBED (MBCOUT , PBCOUT , PHIOUTl , TBCOUT , MDBOUT3 , PDB0UT3 , PHI 
OUT 3 , TDBOUT3 , CPAIR , DT , MS BADS , MH20 , MH20R , MDB 
TANK , 3 , TIME , ADSORB2 ) 

CALL SORPUMP ( PBCOUT , MSPOUT , PSPOUT , TSPOUT , CVAIR , CVBE 

D , CVC02 , DT , HC02 , HEAT , MC02 , MS BADS , MAIR , MSBTA 
NK, PBED, PVC02 , RAIR, RC02 , SPEED, TBED, TIME , VBE 
D, VPUMP, C02DES0RB, MC02P) 

IF ( VALVE. EQ.l) THEN 
MSBOUT = MSBOUT 
MSPOUT = 0 
ENDIF 

CALL DESSBED (MSBOUT, PSBOUT, PHIOUTl, TSBOUT, MDBOUT2 , PDBOUT 
2 , PHIOUT2 , TDBOUT2 , CPAIR, DT , MDBADS , MH20 , MH20 
R , MDBTANK , 2 , TIME , DESORB ) 

ELSE 
HEAT=0 
MCO2R=0 
MC02P=0 
MSPOUT=0 
PSPOUT=0 
TSPOUT=0 
ENDIF 

CALL C02TANK (MSPOUT , PSPOUT , TSPOUT , MTOUT , PTOUT , TTOUT , MTAN 
K , PTANK , TTANK , CVC02 , DT , RC02 , TIME , VTANK) 


RETURN 

END 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c. 


... PROGRAM NAME: DESSBED. FOR 

CALLED FROM: MAIN PROGRAM C02REM. FOR 

... GROUP NAME: CLASSICAL CONTROL 

... CREATED BY: TIM SPRECKER 

DAN WALDECK 

PAUL M. SNIDER 

FEBRUARY 18, 1992 

. .. REVISED: 

. .. CALLS: SUBROUTINE G1Z2 

. .. CALLED FROM: SUBROUTINE DESSBED. FOR 

. .. OUTPUT FILE: NONE 

LIST OF VARIABLES 


... CPAIR : SPECIFIC HEAT OF AIR (kJ/kg*K) 

... DT: DELTA TIME STEP (s) 

... HDBIN : ENTHALPY OF CONDENSATION (kJ/kg) 

... DBLOAD : LOAD ON THE DESSICANT BEDS 

... MDBADS (2,2) : MASS OF WATER ADSORBED BY DESSICANT (kg) 

... MDBIN : MASS OF CABIN AIR INTO THE DESSICANT BEDS (kg/s) 

... MDBOUT : MASS OF CABIN AIR OUT OF THE DESSICANT BEDS (kg/s) 
... MH20: MASS OF WATER IN THE INCOMING CABIN AIR (kg) 

... MH20R: MASS OF WATER REMOVED BY THE DESSICANT (kg) 

MDBTANK ( 2 ) : MASS OF DESSICANT MATERIAL (kg) 

. . . M: INTEGER BED SWITCHING VARIABLE 
... N: INTEGER BED SWITCHING VARIABLE 

... PDBIN : PRESSURE OF CABIN AIR INTO THE DESSICANT BEDS (kPa) 
... PDBOUT : PRESSURE OF CABIN AIR OUT OF THE DESSICANT BEDS 
TO THE BLOWER/ PRECOOLER (kPa) 

PHI IN : RELATIVE HUMIDITY OF CABIN AIR INTO THE DESSICANT 

BEDS (%/100) 

PHI OUT : RELATIVE HUMIDITY OF CABIN AIR OUT OF THE 

DESSICANT BEDS TO THE BLOWER/ PRECOOLER (%/100) 

. .. PVH20: PARTIAL PRESSURE OF WATER VAPOR (kPa) 

. .. Tl: TEMPORARY TEMPERATURE VARIABLE FOR ITERATION (deg K) 

T2 : TEMPORATY TEMPERATURE VARIABLE FOR ITERATION (deg K) 

TDBIN : TEMPEATUTE OF CABIN AIR INTO THE DESSICANT BEDS (K) 

. .. TDBOUT : TEMPERATURE OF CABIN AIR OUT OF THE DESSICANT BEDS 

TO THE BLOWER/ PRECOOLER (K) 

. .. TIME: PRESENT SIMULATION TIME (s) 

. . . W: ABSOLUTE HUMIDITY OF CABIN AIR 


C$NOEXT 

C$NOWARN 


SUBROUTINE DESSBED (MDBIN , PDBIN , PHIIN, TDBIN, MDBOUT , PDBOUT, PH 
& I OUT , TDBOUT , CPAIR , DT , MDBADS , MH20 , MH20R , M 
& DBTANK,N, TIME, DBLOAD) 

C INITIALIZE PARAMETERS 


94 


*>*!*»*> *»*»*>*> 


REAL MDBADS (2,2), MDBTANK ( 2 ) 

DOUBLE PRECISION MH20 , MH20R , MDBIN , MDBOUT 
INTEGER DT , TIME 
M=MOD ( TIME/ 18 0 0,2) 

MDBIN=MDBIN * DT 

IF(N.EQ. 1)THEN 

C CALCULATE PARTIAL PRESSURE OF WATER VAPOR 

PVH20=PHIIN*PSAT (TDBIN-273 ) 

C . . . CALCULATE OMEGA IN KG VAPOR/KG DRY AIR 
W=( . 622*PVH20)/ (PDBIN-PVH20) 

C BREAKDOWN AIR COMPONENTS 

MH20=W*MDBIN/ ( 1+W) 

C ... CALL ZEOLITE SUBROUTINE FOR ADSORPTION AND DESORPTION 
I WEIGHT = 1 

CALL G1Z2 (MDBIN , PDBIN , PHIIN , TDBIN , MDBOUT , PDBOUT , PHIOUT , TD 
BOUT , CPAIR , DT , MDBADS (M+l , 1 ) , MH20 , MH20R , MDBTANK ( 

1) , 1, TIME, DBLOAD, IWEIGHT) 

CALL G1Z2 (MDBOUT , PDBOUT , PHIOUT , TDBOUT , MDBOUT , PDBOUT , PHIOUT , TD 
BOUT, CPAIR, DT, MDBADS (M+l, 2) ,MH20,MH20R, MDBTANK ( 

2 ) , 2 , TIME , DBLOAD , IWEIGHT) 

ELSE IF (N . EQ . 2 ) THEN 
IWEIGHT = 1 
TDBIN=363 

CALL G1Z2 (MDBIN , PDBIN , PHIIN , TDBIN , MDBOUT, PDBOUT , PHIOUT , TD 
BOUT, CPAIR, DT, MDBADS (2 -M, 2) ,MH20,MH20R, MDBTANK ( 

2) ,2, TIME, DBLOAD, IWEIGHT) 

CALL G1Z2 (MDBOUT , PDBOUT , PHIOUT , TDBOUT , MDBOUT , PDBOUT , PHIOUT , TD 
BOUT , CPAIR , DT , MDBADS ( 2 -M , 1 ) , MH20 , MH2 OR , MDBTANK ( 

1 ) , 1 , TIME , DBLOAD , IWEIGHT) 

ELSE 

IWEIGHT = 0 

CALL G1Z 2 (MDBIN , PDBIN , PHIIN , TDBIN , MDBOUT , PDBOUT , PHIOUT , TD 
BOUT , CPAIR , DT , MDBADS (M+ 1 , 2 ) , MH20 , MH20R , MDBTANK ( 

2) ,2, TIME, DBLOAD, I WEIGHT) 

c write (*,*)mh2or,phiin,phiout 

ENDIF 

MDBIN=MDBIN/DT 
MDBOUT=MDBOUT/ DT 

RETURN 

END 
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SUBROUTINE G1Z2 ( MDBIN , PDBIN , PHI IN , TDBIN , MDBOUT , PDBOUT , PHIOU 
f T , TDBOUT , CPAIR , DT , MDBADS , MH20 , MH20R , MDBTANK 

& , N , TIME , DBLOAD , IWEIGHT) 

C INITIALIZE PARAMETERS 

REAL MDBADS , MDBTANK 

DOUBLE PRECISION MH20,MH20R, MDBIN, MDBOUT 
INTEGER DT, TIME 

C . . . DETERMINE CURRENT TANK LOAD 
ICOUNT=0 

DBLOAD=MDBADS /MDBTANK 

C . . . FIGURE AIR CHARACTERISTICS AT EQUILIBRIUM WITH DESICCANT 

IF(N.EQ.l) PHIOUT=DBLOAD/ .5263 
IF (N.EQ. 2 .AND. DBLOAD. LE. . 17) PHIOUT= . 4 *DBLOAD 
IF(N.EQ. 2. AND. DBLOAD. GT. .17) PHIOUT= . 068+40* (DBLOAD- 17) 
T1=TDBIN ; 

C ... CALCULATE PARTIAL PRESSURE OF WATER VAPOR 
10 PVH20=PHI0UT*PSAT (Tl-273 ) 

W=( . 622*PVH20)/ (PDBIN-PVH20) 

C ... CALCULATE WATER REMOVED 

IF ( IWEIGHT. EQ. 1) THEN 

MH20R=MH20-W* (MDBIN-MH20) 

ELSE 

MH20R=5* (MH20-W* (MDBIN-MH20) ) 

ENDIF 


C . 


C . 
C . 


C . 


ADD ENERGY OF EVAPORATION TO AIR AND FIND 
HDBIN=MH20R* (2502-2 . 389* (TDBIN-273) ) 
T2=TDBIN+HDBIN/ (MDBIN-MH20R) /CPAIR 


NEW TEMPERATURE 


ITERATE TO FIND EQUILIBRIUM TEMPERATURE 
IF (ABS (T2-T1) . GT . 1 ) THEN 
. WEIGHT THE AVERAGE TOWARD T1 FOR STABILITY 
Tl=(3*Tl+T2)/4 
ICOUNT=ICOUNT+l 

QUIT WITH CURRENT T1 AFTER 50TH TRY 
IF(ICOUNT. LT. 50) GOTO 10 
ENDIF 


TDBOUT=Tl 

PDBOUT=PDBIN 

C . . . ADD REMOVED WATER TO BED 
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MDBADS =MDBADS +MH 2 OR 


C FIGURE NEW MAKEUP OF AIR 

MH20=MH20-MH20R 
MDBOUT=MDBIN-MH2 OR 

RETURN 

END 

C WATER VAPOR SATURATION PRESSURE (kPa) AT TEMPERATURE (deg C) 

FUNCTION PS AT (T) 

PSAT= . 3972+ . 0629*T+ . 00 109 9 *T** 2+ . 00001705*T**3+ . 0000006192*T**4 
END 
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PROGRAM NAME: BLOWCOOLFOR 

c ::: SESF-ST* <SSiSS B SS»S a »“-'« 

C ... CREATED BY: CARL ALBRECHT 

c ROGER BURJES 

c PAUL M. SNIDER 

^ * MARCH 12, 1992 

C ... REVISED: 

C . . . SUBROUTINE CALLS : NONE 

C ... OUTPUT FILE: NONE 

C * LIST 0F VARIABLES 

C ... CPH20: SPECIFIC HEAT* OF* AIR* 

c ::: (k3) 

C ... PBCIN : PRESSURE OF DESSICANT BED aIp T m™ 

C . . . PBCOUT : PRESSURE OF AIROOTOPto^LS THE BL0WER (kPa) 

c QDOT: HEAT TRANSFER rrtuppm 2f^ r THE PREC0 °LER (kPa) 

c . . . TBCIN: TEMPERATURE OF^SSICANT PREC00LER (W/s) 

c ::: gggjggs £ S“«* 7 ™ E (K) 

§ : ! EL0C y TY , 0 F ,“ bin wwS’iS (K) 



C$NOEXT 

C$NOWARN 

SUBROUTINE BLOWCOOL(MBCIN, PBCIN, TBCIN, MBCOUT, PBCOUT TBCOUT 

OPH20, EPSILON, PAIR, TWAIE^ VEMR^ ' 

DOUBLE PRECISION MBCIN, MBCOUT ~~ 

C CALCT “TE HEAT TRANSFER BETWEEN BLOWER AIR AND COOLANT 
mbcot^mb“n* CPH2 °* MBCIN * (tbcin -TWATER) 

PBCOUT = PBCIN 


C . . . FIND NEW BLOWER AIR TEMPERATURE 
TBCOUT=TBCIN-QDOT/ ( CPH20*MBC0UT) 

RETURN 

END 
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C$NOEXT 

C$NOWARN 

C$INCLUDE FUNCZ5A. FOR 


SUBROUTINE SORBED (MSBIN, PSBIN, TSBIN,MSBOUT, PSBOUT, TSBOUT, CP 
& AIR , CVBED , DT , HC02 , MS BADS , MSBTANK , PVC 

& 02 IN , PVC020UT , TBED , TIME , SBLOAD , MC02R) 

REAL MSBTANK (3) ,MSBADS(2,3) ,TBED(2,3) 

DOUBLE PRECISION MC02 ,MC02R, MSBIN, MSBOUT 
INTEGER DT, TIME 
M=MOD ( TIME/ 1800,2) 

MSBIN=MSBIN*DT 

MC02=1 . 519*PVC02IN/PSBIN*MSBIN 

C . . . DETERMINE CURRENT TANK LOAD 

S BLOAD= 1 0 0 * MS BADS (M+l , 3 ) /MSBTANK ( 3 ) 

C • • • FIGURE AIR CHARACTERISTICS AT EQUILIBRIUM WITH SORBENT 
PVC020UT=Z 5 A (SBLOAD , TBED (M+ 1 , 3 ) -2 7 3 ) 

W=1 . 519*PVC020UT/ ( PSBIN-PVC020UT) 

C CALCULATE C02 REMOVED 

MC02R=MC02-W* (MSBIN-MC02 ) 

IF (— MC02R. GT . MS BADS (M+l , 3 ) ) MC02R=-MSBADS (M+l , 3 ) 

C . . . ADD REMOVED C02 TO BED 

MS BADS (M+l , 3 ) =MSBADS (M+l , 3 ) +MC02R 

C . . . FIGURE NEW MAKEUP OF AIR 

MS BOUT=MS BIN-MC02R 

C ... ADD ENERGY OF EVAPORATION TO AIR+BED AND FIND NEW 
C TEMPERATURE 


TSBOUT= (MC02R*HC02+MSB0UT*CPAIR*TSBIN+ (MSBTANK ( 3 ) +MSBA 
& DS (M+l , 3 ) ) *CVBED*TBED (M+l , 3 ) ) / (MSBOUT*CPAIR+ (MSBTANK 

& (3) +MSBADS (M+l, 3) ) *CVBED) 

TBED ( M+ 1 , 3 ) =TS BOUT 
PSBOUT=PSBIN 
MSBIN=MSBIN/DT 
MSBOUT=MS BOUT / DT 


RETURN 

END 
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C$NOEXT 

C$NOWARN 


SUBROUTINE SORPUMP ( PSPIN , MSPOUT , PSPOUT , TSPOUT , CVAIR , 

& CVBED , CVC02 , DT , HC02 , HEAT , MC02 , MS BADS , MAIR , 

& MSBTANK, PBED , PVC02 , RAIR, RC02 , SPEED , TBED , 

& TIME, VBED, VPUMP, SBL0AD,MC02P) 

REAL MSBTANK ( 3 ) , MS BADS (2,3), TBED (2,3) 

DOUBLE PRECISION MSPOUT, MC02P, MAIR, MC02 ,MC02R 
INTEGER DT , TIME 
M=MOD(TIME/1800, 2) 

C . . . DETERMINE CURRENT TANK LOAD 

S BLOAD= 1 0 0 *MS BADS ( 2 -M , 3 ) /MS BTANK ( 3 ) 

IF (M. NE . MOD ( (TIME-DT) /1800 , 2 ) ) THEN 
PBED=PSPIN 
HEAT=12 . *DT 

C ... FIGURE INITIAL AIR CHARACTERISTICS AT EQUILIBRIUM WITH 
C . . . SORBENT 

PVC02=Z5A(SBL0AD,TBED(2-M, 3) -273) 

MC02=PVC02 * VBED/RC02 /TBED ( 2 -M , 3 ) 

MAIR= (PBED-PVC02 ) *VBED/RAIR/TBED (2-M, 3) 

ENDIF 

IF (TBED (2-M, 3) ,GT. 478 .AND. HEAT. GT. 0 . ) HEAT=0.*DT 

IF (TBED (2-M, 3) .GT. 368 .AND. MOD (TIME, 1800) .GE.1480) HEAT=-12 . * 

& DT 


C ... ADD ENERGY INPUT TO AIR+BED AND FIND NEW TEMPERATURE AND 
C . . . PVC02 


TBED (2-M, 3 ) =TBED (2-M, 3 ) +HEAT/ (MAIR*CVAIR+MC02*CVC02+ (MSBTAN 
& K(3) +MSBADS (2-M, 3) ) *CVBED) 

PVC02=Z5A (SBLOAD , TBED ( 2 -M , 3 ) -2 7 3 ) 

C ... CALCULATE C02 REMOVED FROM BED 

MC02R=PVC02 *VBED/RC02 /TBED ( 2 -M , 3 ) -MC02 

C SUBTRACT REMOVED C02 FROM BED 

MS BADS ( 2 -M , 3 ) =MS BADS ( 2 -M , 3 ) -MC02R 
MC02 =MCO 2 +MCO 2 R 

C ADD ENERGY OF EVAPORATION TO AIR+BED AND FIND NEW 

C . . . TEMPERATURE 


TBED (2-M, 3 ) =TBED ( 2-M, 3 ) -MC02R*HC02/ (MAIR*CVAIR+MC02*CVC02+ ( 
& MSBTANK(3) +MSBADS (2-M, 3) ) *CVBED) 
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C RUN CO 2 PUMP 

MC02P=MC02*(VPUMP*SPEED*DT/60/(VBED+VPUMP*SPEED*DT/60 ) ) 

MSPOUT= (MC02+MAIR) * (VPUMP*SPEED*DT/60/ (VBED+VPUMP*SPEED* 

& DT*60 . ) ) 

PSP0UT=PVC02+MAIR*RAIR*TBED(2-M, 3) /VBED 
TSPOUT=TBED ( 2 -M , 3 ) 

MC02=MC02* ( 1-VPUMP*SPEED*DT/ 60/ (VBED+VPUMP*SPEED*DT/60 . ) ) 
MAIR=MAIR* (l-VPUMP*SPEED*DT/60/ (VBED+VPUMP*SPEED*DT/60 . ) ) 

C . . . DETERMINE CURRENT TANK LOAD 

SBLOAD=100*MSBADS (2-M, 3 ) /MSBTANK ( 3 ) 

C . . . FIGURE AIR CHARACTERISTICS AT EQUILIBRIUM WITH SORBENT 

PVC02=Z5A (SBLOAD, TBED (2-M, 3 ) -273 ) 

C CALCULATE C02 REMOVED FROM BED 

MC02R=PVC02*VBED/RC02/TBED (2-M, 3 ) -MC02 

C ... SUBTRACT REMOVED C02 FROM BED 

MS BADS (2-M, 3 ) =MSBADS (2-M, 3 ) -MC02R 
MC02=MC02+MC02R 

. . . ADD ENERGY OF EVAPORATION TO AIR+BED AND FIND NEW 
. . . TEMPERATURE 

TBED ( 2-M , 3 ) =TBED ( 2 -M , 3 ) -MC02R*HC02/ (MAIR*CVAIR+MC02 *CVC02+ ( 
& MSBTANK(3) +MSBADS (2-M, 3) ) *CVBED) 

PBED=PVC02+MAIR*RAIR*TBED (2-M, 3 ) /VBED 

RETURN 
END 
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c 

• • • 

• m m 

PROGRAM NAME 

: C02TANK. FOR 

c 

• • • 

CALLED FROM: 

MAIN PROGRAM C02REM. FOR 

c 

• • • 

GROUP NAME: 

CLASSICAL CONTROLS 

c 

• ♦ • 

CREATED BY: 

DANIEL T. WALDECK 

c 



PAUL M. SNIDER 
MARCH 29, 1992 

c 



c 

• # # 

REVISED: 

c 

• • • 

CALLS : NONE 


c 

c 

• • • 

t • • I 

OUTPUT FILE: 

NONE 


C$NOEXT 

C$NOWARN 


SUBROUTINE C02TANK (MTIN , PTIN , TTIN , MTOUT , PTOUT , TTOUT , MTANK , P 
& TANK , TTANK , CVC02 , DT , RC02 , TIME , VTANK) 

DOUBLE PRECISION MTIN , MTOUT , MTANK 
INTEGER DT , TIME 

C . . . CALCULATE CURRENT TANK CONDITIONS 

TTANK=TTANK+MT IN* CVCO 2 *TTIN/ (MTIN+MTANK) 

PTANK=MTANK*RC02 *TTANK/VTANK 

C . . . CALCULATE NEW MASS OF THE TANK 

MTANK=MTANK+MT I N -MTOUT 

C . . . CALCULATE TANK OUTPUT CONDITIONS 

TTOUT=TTANK-MTOUT * CVCO 2 *TTANK/ ( MTANK ) 

PTOUT= PT AN K 

RETURN 

END 
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5.8 Graphical Interface Code 


/* 

♦Filename : show . h 
* 

* Creator: Carl Albrect 

* 

* Created: March 31, 1992 

* Modified: May 9, 1992 - Robert Swenson - added bed switch 

* — 

* Description: 

* types and constants for the show program 


struct diag_data { 
int C02 , 

Tank_Pres , 

Pump_Speed , 

Absorb, 

Desorb, 

C02_Ab, 

C02_De , 

sys, 

fan; 

int oldC02 , /* what the graph was last time step */ 

oldTank_Pres , 
oldPump_Speed , 
oldAbsorb, 
oldDesorb, 
oldC02_Ab, 
oldC02_De ; 

int bed_switch, /* o if top bed = desorb, 1 otherwise */ 

old bed ; 
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/*\ 

* Filename: show.c 

* 

* Created: March 15, 1992 - converted from bars.c 

* Modified: May 9, 1992 Robert Swenson - added code . bed_switch 

* 

* creator: Robert Swenson 

* 

* Description: 

* create a window, list NASA/USRA ADT members, show overview of 

* carbon dioxide removal system, and show pressure bars/graphs 

W 


# include <X11/Xlib.h> 

# include <X11/Xutil .h> 

# include <Xll/Xos.h> 

#include <Xll/Xatom.h> 

# include "show.h" 

# include <stdio.h> 

#include "icon_bitmap" 

#def ine BITMAPDEPTH 1 
#define TOO_SMALL 0 
#define BIG_ENOUGH 1 

Display *display; 
int screen_num; 

static char *progname; /* name this program was invoked by (argv[0]) */ 

void main(argc, argv) 
int argc; 
char **argv; 

{ 

Window win; 

unsigned int width, height; /* window size */ 

int x, y; /* window position */ 

unsigned int border_width =4; /* four pixels */ 

unsigned int display_width, display_height ; 

unsigned int icon_width, icon_height; 

char *window_name = ''Simulation Demonstration"; 

char *icon_name = "Simulation Demo"; 

Pixmap icon_pixmap ; 

XSizeHints size_hints; 

XlconSize *size_list; 
int count; 

FILE *my_data; 

XEvent report; 

GC gc; 


104 



XFontStruct *font_info; 
char *display_name = NULL; 

int window_size = BIG_ENOUGH; /* or TOO_SMALL to display contents */ 
struct diag_data code; 

progname = argv[0]; 

/* connect to X server */ 

if ( (display=XOpenDisplay (display_name) ) == NULL ) 

{ 

(void) fprintf( stderr, ”%s: cannot connect to X server %s\n", 
progname, XDisplayName (display_name) ) ; 
exit( -1 ) ; 

) 

/* get screen size from display structure macro */ 
screen_num = DefaultScreen (display) ; 
display_width = DisplayWidth (display , screen_num) ; 
display_height = DisplayHeight (display , screen_num) ; 

/* Note that in a real application, x and y would default to 0 
* but would be settable from the command line or resource database. 
*/ 

x = y = 0; 

/* size window with enough room for text */ 
width = 640, height = 460; 

/* create opaque window */ 

win = XCreateSimpleWindow (display , RootWindow (display, screen_num) , 
x, y, width, height, border_width, BlackPixel (display , 
screen_num) , WhitePixel (display, screen_num) ) ; 

/* Get available icon sizes from Window manager */ 

if (XGetlconSizes (display , RootWindow (display, screen_num) , 
&size_list, &count) == 0) 

(void) fprintf ( stderr, 

"%s: Window manager didn't set icon sizes - using default. \n", 
progname) ; 

else { ; /* A real application would search through size_list 

* here to find an acceptable icon size, and then 

* create a pixmap of that size. This requires 

* that the application have data for several sizes 

* of icons. */ 

) 

/* Create pixmap of depth 1 (bitmap) for icon */ 

icon_pixmap = XCreateBitmapFromData (display , win, icon_bitmap_bits, 
icon_bitmap_width, icon_bitmap_height) ; 
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size_hints. flags = PPosition J PSize j PMinSize; 
size_hints.min_width = width; 
size_hints.min_height = height; 

{ 

XWMHints wm_hints; 

XClassHint class_hints; 

/* format of the window name and icon name args has changed in R4 */ 
XTextProperty windowName, iconName; 

if (XStringListToTextProperty (&window_name, 1, SwindowName) == 0) { 
(void) fprintf ( stderr, 

"%s: structure allocation for windowName failed. \n", progname); 
exit(-l) ; 

} 

if (XStringListToTextProperty (&icon_name, 1 , &iconName) == 0 ) { 

(void) fprintf ( stderr, 

"%s: structure allocation for iconName failed. \n", progname); 
exit(-l) ; 

> 

wm_hints . initial_state = NormalState; 

wm_hints . input = True ; 

wm_hints. icon_pixmap = icon pixmap; 

wm_hints. flags = StateHint j IconPixmapHint | InputHint; 

class_hints . res_name = progname ; 
class_hints.res_class = ''Simulation Demo"; 

XSetWMProperties (display, win, SwindowName, &iconName, 
argv, argc, &size_hints, &wm_hints, 

&class_hints) ; 

) 

/* Select event types wanted */ 

XSelectlnput (display , win, ExposureMask j KeyPressMask ] 
ButtonPressMask ] StructureNotifyMask) ; 

load_font ( &font_info) ; 

/* create GC for text and drawing */ 
getGC(win, &gc, font_info) ; 

/* Display window */ 

XMapWindow( display, win) ; 

/* open file */ 

if ( (my_data = fopen ("data. sim" , "r") ) == NULL) 

fprintf (stderr, "bar: unable to open data.simXn") ; 
exit (-1) ; 
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} 


/* init the structure */ 
code.C02 = 0; 
code . Tank_Pres = 0; 
code . Pump_Speed = 0; 
code. Absorb = 0; 
code . Desorb = 0 ; 
code.C02_Ab = 0; 
code . C02_De = 0 ; 
code. fan = 0; 
code . bed_switch = 0; 

/* init the window */ 

draw_text (win, gc, font_info, width, height, code.bed_switch) ; 
draw_graphics (win, gc, width, height, &code) ; 


/* get events, use first to display text and graphics 
/* To run through the data file infinity times */ 
while (1 == 1) 

{ 

if (feof (my_data) ) 

{ 

fseek(my_data, 0 , 0) ; 

} 


/* don't forget the old data */ 
code.oldC02 = code.C02; 
code . oldTank_Pres = code . Tank_Pres ; 
code . oldPump_Speed = code . Pump_Speed ; 
code. oldAbsorb = code. Absorb; 
code.oldDesorb = code. Desorb; 
code.oldC02_Ab = code.C02_Ab; 
code. oldC02_De = code.C02 De; 


*/ 


if (! XPending (display) ) /* if nothing to do then just draw */ 

{ 

/* get the next piece of data */ 
fscanf (my_data, "%i\t%i\t%i\t%i\t%i\t%i\t%i\t%i\n" , &code.C02 , 
&code.Tank_Pres, & c ode. Pump_S peed, &code. Absorb, 

&code. Desorb, Scode . C02_Ab, &code . C02_De , Scode.sys); 

code . bed_switch = 0 ; 
if (code.bed_switch != code. old_bed) 

{ 

/* erase the window and re-draw with the new beds */ 

) 


draw_data( win, gc, width, height, &code) ; 

) 

else 

{ /* else update the win or whatever */ 

XNextEvent (display , Sreport) ; 
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switch (report. type) { 
case Expose: 

/* unless this is the last contiguous expose, 

* don't draw the window */ 
if (report . xexpose . count != 0) 

break; 

/* if window too small to use */ 

/*if (window_size == T00_SMALL) 

TooSmall (win, gc, font_info) ; 
else*/ { 

/* place text in window */ 

draw_text (win, gc, font_info, width, height, 
code . bed_switch) ; 

/* place graphics in window, */ 
draw_graphics (win, gc, width, height, &code) 

break ; 

case Conf igureNotify : 

/* window has been resized, change width and 

* height to send to draw_text and draw_graphics 

* in next Expose */ 

width = report. xconfigure. width; 
height = report . xconfigure. height ; 
if ((width < size_hints .min_width) Jj 

(height < size_hints.min_height) ) 
window_size = T00_SMALL; 

else 

window_size = BIG_ENOUGH; 
break; 

case ButtonPress: 

/* trickle down into Keypress (no break) */ 
case KeyPress: 

XUnloadFont (display , font_info->f id) ; 

XFreeGC (display , gc) ; 

XCloseDisplay (display) ; 
exit (1) ; 
default: 

/* all events selected by StructureNotifyMask 

* except Conf igureNotify are thrown away here, 

* since nothing is done with them */ 
break; 

) /* end switch */ 

) /* XPending */ 

) /* end while */ 

f close (my _data) ; /* we're done with the data */ 

char ch; 

ch = getchar ( ) ; 


} 
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} /* main */ 


getGC(win, gc, font_info) 

Window win; 

GC *gc ; 

XFontStruct *font_info; 

{ 

unsigned long valuemask = 0; /* ignore XGCvalues and use defaults */ 

XGCValues values; 

unsigned int line_width = 3; 

int line_style = LineSolid; 

int cap_style = CapRound; 

int join_style = JoinRound; 

int dash_offset = 0; 

static char dash_list[] = {12, 24}; 

int list_length = 2; 

/* Create default Graphics Context */ 

*gc = XCreateGC( display, win, valuemask, Svalues) ; 

/* specify font */ 

XSetFont (display, *gc, font_info->f id) ; 

/* specify black foreground since default window background is 
* white and default foreground is undefined. */ 

XSetForeground (display, *gc, BlackPixel (display , screen_num) ) ; 

/* set line attributes */ 

XSetLineAttributes (display, *gc, line_width, line_style, 
cap_style, join_style) ; 

} /* getGC */ 


load_f ont ( f ont_info) 

XFontStruct **font_info; 

{ 

char *fontname = ''SxlS” ; 

/* Load font and get font information structure. */ 
if ((*font_info = XLoadQueryFont (display , fontname) ) == NULL) 

{ 

(void) fprintf (stderr , "%s: Cannot open 9x15 font\n",progname) ; 
exit( -1 ) ; 

> 

} /* load_font */ 


TooSmall (win, gc, font_info) 
Window win; 
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GC gc; 

XFontStruct *font_info; 

{ 

char *stringl = "Too Small"; 
int y_offset, x_offset; 

y_of f set = font_info->ascent + 2 ? 
x_offset = 2; 

/* output text, centered on each line */ 
XDrawString (display , win, gc, x_offset, y_offset, 
strlen (stringl) ) ; 

} /* TooSmall */ 


stringl , 
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/*\ 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


w 


Filename: draw.c 

Created: March 15, 1992 - converted from bars.c 


creator: Carl Albrecht 

Modified: Robert Swenson - April 
Modified: RAS - May 6 - reduced 


3, 1992 - added draw_data 

window/drawing size to fit on PC 


Description: 

Draw in the window 


# include <X11/Xlib.h> 

# include <X11/Xutil .h> 

# include <Xll/Xos.h> 

# include <Xll/Xatom.h> 

# include <stdio.h> 

# include "show.h" 

extern Display *display; 


#def ine 
#def ine 
*/ 

#define 
#define 
#def ine 
#define 
#define 
#define 
#define 
#define 


Nnmstrinas 30 /* number of strings to print on the screen */ 

?extBa?Width 280 /* width in pixels o t region containing text 


TextTab 9 

VOffset 10 

StatBox 225 

StatBarTitles 420 

NumSegments 70 /* # 

WordsInPic 25 


/* # of pixels to tab over for uncentered */ 
/* move down x pixels before typing */ 

/* y coord of the Cabin Status box */ 

/* "base” of status bars */ 
of unconnected lines to draw */ 


/ft ^ OI UIIUUIUIC^ ( 

_ __ /* | of words to print into the picture / 

Ba num a Ws XH 4 19ht /* num°of names to switch when beds switch */ 


draw_text (win , gc, font_info, win_width, win_height, bed_flag) 

Window win; 

GC gc; 

XFontStruct *font_info; 

unsigned int win width, win_height, 

int bed_flag; 

{ static char *strings [NumStrings] = { 

"NASA/USRA” , 

"ADT" , 

"at", 

"Kansas State University", 

"Carl Albrecht Roger Burjes" , 

"Dan Waldeck Stan Kroeker", 
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"Maury Wilmoth 
"Mike Honas”, 
"Mike Brockway 
"Paul Snider 
"Robert Swenson 


Tim Sprecker” , 

Dr. Cogley", 

Dr . Gustafson” , 
Robert Young” , 


"Cabin Status Panel”, 
"Carbon Tank Fan", 

"Dioxide Press Speed”, 


"Fan” , 
"Dryer” , 
"Pump" , 

II II 

t 

II »l 

"Desorb" , 
"Tank" , 
"Pre" , 
"Cooler" , 
"CO" , 

" 2 ", 


/*14*/ 

/*16*/ 

/*18*/ 

/* 20 */ 

/* 22 */ 


/* words in picture */ 


/*24*/ 

"Rehumidifier" , 

"Absorb", /*26*/ 

"Storage" 

static struct ^ to *' 

/* strings [] indexes and just run through / 

/* the array printing the words*/ 
int x /* coord to put the word at */ 


the_word; /* which strings [] to put at x,y */ 

} the coords [WordsInPic] = { 

7410,320,14}, 

{560,100,16}, 

{545,240,23}, 

{565,245,24} , 

{545,355,23} , 

{565,360,24} , 

{465,145,20} , 

{470,335,21} , 

{460,350,22} , 

{470,105,23} , 

{490,110,24} , 

{455,125,27} 

bed_coords[ 2 ][num_beds] = { { */ 

{320,240,15} , 

{545,260,19} , 

{545,375,26} , 

{295,450,25} 


/* beas switcn . , 


/* Desorb */ 
/* Adsorb */ 
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>» 

{ 

{320,240,25} , 

{545,260,26} , 

{545,375,19} , 

{295,450,15} 

} 

}; 

int i; 
int len; 

int fontjaeight = font_info->ascent + f ont_inf o->descent ; 


*/ 


for ( i= 0 ; i<4 ; i++) /* CENTER THE TITLES */ 

efrion^trinasril ) ; /* len for XTextWidth and XDrawString */ 

i!Sth=X^ width for centering 

XDrawString (display, win, go, (TextBarWidth - width] 1/2, 

VOffset + (i + 1) * (1 + font_height) , stnngs[i], len; , 

} 

for ( i=4 ; i<ll ; i++) /* SHOW THE NAMES */ 

'len = strlen ( strings [ i ] ) ; /* len for XTextWidth and XDrawString */ 

XDrawString (display, win, gc, ^extTab, na sri1 len) ; 

VOffset + (i + 2) * (1 + font_height) , stnngs[i], ienj 

} 

len = str len ( strings^ 11 ] ^ ATUS / len ‘for XTextWidth and XDrawString V 

XDraSs«!^g"dtsplf?? t Sin!°4S? r (?lrt MrWiSth' - width) /2 , 

StatBox + 2 * VOffset, stnngs[ll], len), 
for (i=12 ; i<14 ; i++) /* SHOW STATUS BAR TITLES */ 

{ ien = strlen ( strings [ i ] ) ; /* len for XTextWidth and XDrawString */ 

XDrawString (display, win, gc, ^extTab, . ■ hn 

StatBarTitles + (1 - 11) * (1 + font_heigh ), 

strings [i], len); 

for (i=0;i<WordsInPic;i++) /* PUT WORDS IN THE PICTURE */ 

ion — efri fin ( str incjs r the coords [ 1 } • the_word ] ) , _ . 

width= XTextWidth ( f ont_inf o , strings [ the coords [ i ] .the_word] , len) . 
XDrawString (display , win, gc, the_coords [ 1 ] . x, len) * 

the_coords [ i ] . y , strings [ the_coords [ 1 ] . the_word ] , ) , 

} 


for ( i=0 ; i<num_beds ; i++) /* PUT IN THE bed names */ 
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^len = strlen ( strings [bed_coords[bed_f lag] [i] .the_word] ) ; 
width=XTextWidth(font_info, strings [bed_coords[bed_f lag] [ 1 ] .the_word] , 

len) ; 

/* INCASE the beds just switched, clear the area where this */ 
/* word is about to go to avoid overwrite */ 

/* XClearArea ( ) ; */ 


XDrawString (display , win, gc, bed_coords[bed_flag] [ 1 ] .x, 
bed_coords[bed_flag] [i] .y, 

strings [bed_coords [bed_f lag] [i] .the_word] ,len ) , 

} 

} /* draw_text */ 


#define 
#define 
#def ine 
#define 
#def ine 
#def ine 
#define 
#define 
#def ine 
#define 
#define 
#define 
#define 
#define 
#define 


/* x coord of C02 tank bar graph */ 

/* x coord of CO 2 tank pres bar graph */ 

195 /* x coord of pump rpm bar graph */ 

/* x coord of Deisecant bed box in diagram */ 

x~fan 411 /* x coord of the fan */ 
xjre 460 /* X coord of precooler */ 

x C02des 550 /* x coord of Zeolite beds */ 

BarWidth 10 /* width of status bars */ 

BoxWidth 47 /* width of the diagrams boxes */ 

C Width 28 /* width of the fan circle */ 

C _ Height 28 /* height of the fan circle */ 

Boxheight 53 /* height of boxes in pixels */ 

y_box 2 
y_rpm 0 
y_pres 1 


x_C02 60 

x_tank 120 
x_rpm 


draw_graphics (win, gc, window_width, window_height, code) 
Window win; 

GC gc; . ^ . . . 

unsigned int window_width , window_height ; 
struct diag_data *code; 


int x, y, oldy[3]; 
int height; 
int width ; 

FILE *my_data ; 

int i, j, Time_Step, rpm; 

f lost Pros ? 

static XSegment segments [NumSegments] = { 

{ 0, StatBox, TextBarWidth, StatBox}, 

{ TextBarWidth, 0, TextBarWidth, 0), 

{ 0, StatBarTitles, TextBarWidth, StatBarTitles } , 
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{ 280,190,320,190}, 
{ 280,200,320,200}, 
{ 280,400,320,400}, 
{ 280,410,320,410}, 
{ 370,180,550,180}, 
( 370,190,550,190}, 
{ 370,200,400,200}, 
{ 370,210,390,210}, 


/* pipe A top */ 
/* A Bottom */ 

/* B Top */ 

/* B Bot */ 

/* G Top */ 

/* G Bot */ 

/* D Top */ 

/* D Bot */ 


{ 370,390,390,390}, /* C Top */ 

{ 370,400,400,400}, /* C Bot */ 

{ 370,410,550,410}, /* H Top */ 

{ 370,420,550,420}, /* H Bot */ 

{ 390,210,390,390}, /* E left */ 

{ 400,200,400,280}, /* E Right Top */ 
{ 400,400,400,290}, /* E Right Bot */ 
{ 400,280,410,280}, /* F Top */ 

{ 400,290,410,290}, /* F Bot */ 

{ 440,280,460,280}, /* I Top */ 

{ 440,290,460,290}, /* I Bot */ 

{ 510,280,520,280}, /* J Top */ 

{ 510,290,520,290}, /* J Bot */ 

{ 520,280,520,200}, /* K Top Left */ 

{ 520,290,520,400}, /* K Bot Left */ 

{ 520,200,550,200}, /* L Top */ 

{ 520,400,550,400), /* M Bot */ 

{ 530,210,530,390}, /* K Right */ 

{ 530,210,550,210}, /* L Bot */ 

{ 530,390,550,390), /* M Top */ 

{ 600,200,610,200}, /* N Bot */ 

{ 600,190,610,190}, /* N Top */ 

{ 600,410,620,410), /* O Bot */ 

{ 600,400,610,400), /* 0 Top */ 

{ 610,400,610,200}, /* P Bot Left */ 

{ 425,270,425,300}, /* Fan blades */ 
{ 410,285,440,285}, /* Fan Blades */ 
{ 416,275,435,295}, /* Fan Blades */ 
{ 416,295,435,275}, /* Fan Blades */ 
{ 610,190,610,80}, 

{ 620,410,620,70), 

{ 610,80,580,80), 

{ 620,70,590,70}, 

{ 560,60,510,60}, 

{ 570,50,510,50), 

{ 450,60,370,60}, 

{ 450,50,370,50}, 

{ 370,40,370,70), 

{ 300,40,300,70}, 

{ 350,20,320,20}, 

{ 350,90,320,90), 

{ 370,40,350,20}, 

{ 370,70,350,90}, 


/* f Top i-ieiu *•/ 


/* 

P Right */ 

/* 

Q Bot 

*/ 

/* 

Q Top 

*/ 

/* 

R Bot 

*/ 

/* 

R Top 

*/ 

/* 

S Bot 

*/ 

/* 

S Top 

*/ 

/* 

Draws 

the 

/* 

Draws 

the 

/* 

Draws 

the 

/* 

Draws 

the 

/* 

Draws 

the 

/* 

Draws 

the 


octagon/ OGA 
octagon/OGA 
octagon/OGA 
octagon/OGA 
octagon/OGA 
octagon/OGA 


*/ 

*/ 

*/ 

*/ 

*/ 

*/ 
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}; 


< 3oo 0 ;7o 0 :3 3 2 2 o 0 : 9 2 o 0 ’; ;; octagons v 

c SS:n:SS:%!i HE£^r° GA v 

( 480^280^485^270} ,7 *^ raws the Pre?ooler*/ 

{ 485,270,490,280},' 

{ 460,330,480,300), 

{ 480, 300, 485, 310} , 

{ 485,310,490,300}, 

{ 490,300,510,330), 

{ 490,280,510,250} 


can't 


*/ 


segments [ 1 ] . y2 - window.height; windo„_height is . so ^ 

be used in the aggragrate assignment */ 

/* draw the picture */ 

“g Bents(dlsplay< Wln ' ~> /. draw pipes 

XDra^Sanfle s ' 170 . 47 , Boxheight, ; /, aa 

ssssssi:® b e b v 

XDrawRectangle (display,' win,'^' £££“' ™', ^ U* CC V 

/* This draws the fan */ ' 

XDrawArc (display, win, gc, 410, 270,30,30,0,360 * 64) . 

/* This draws the pump */ 

... ... 

vn Th m S draws the C02 Storage Tank */ 

XDrawArc (display, win, gc, 450^25^60^60^0,360 . 64) , 

/* When drawing/redrawing the whole pic, draw complete level bars 

cSIef?oldCO^ t 0 ; he Changed Part */ 

code->oldTank_Pres = 0; 
code->oldPump_Speed = o • 
code->oldAbsorb = 0 ; 
code— >oldDesorb = o ; 
code->oldC02_Ab = 0; 
code->oldC02_De = 0; 

) /* drawlgraphics , , Widtd,he ^ ght ' code ' X) ■' 


*/ 
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draw_data ( win , gc , width , height , code , y ) 

Window win; 

GC gc ; 

int width ; 

int height; 

struct diag_data *code ; 

int y ; 

/* draw the moving bars */ 

if (code->C02 > BarGraphMaxHeight) 

K code->C02 = BarGraphMaxHeight; 

height = BarGraphMaxHeight - ^ od |:^BarTitles - BarGraphMaxHeight, 
XClearArea( display, win, x_C02 , StatBariici 
BarWidth, height, 0) ; 

y = statBarTitles - code->C02 ; 

height = code->C02; BarWidth, height); 

XFillRectangle (display, win, gc, x_co^, y, 


if ( code->Tank_Pres > BarGraphMaxHeight) . 

code->Tank_Pres = BarGraphMaxHeigh , 

XClearArea (display , win, x_tank, StatBarTitles - 

BarGraphMaxHeight^.^ BarGraphMaxHeigh t - code->Tank_Pre S 

y = StatBarTitles - code->Tank_Pres ; 


0 ) ? 


height = code->Tank Pres; v tank V, BarWidth, height); 

XFillRectangle (display, win, gc, x_tank, y, 


if (code->Pump Speed > BarGraphMaxHeight) 

( code->^5_Speed - BarGraphMaxHeight; 


i— x'lTU.lUfce' »_F£/WWe— 

XClearArea (display^win^x_rpm^statBarTitles^-^BarGraphMaxHeight, 

y = StatBarTitles - code->Pump_Speed; 

height = code->Pump Speed; BarWidth, height); 

XFillRectangle (display, win, gc, x_rpm, y, 


y* draw the box fillings */ 

if (code->Absorb > code->oldAbsorb) 

( 


117 



C Sfn-rn^; 2 , y. BoxWidth, height, » 


else if ( code->Absorb < code->oldAbsorb) 


he ight 2 =" code->oldAbsorb ^^““Siwidth- 2 , height, 0) ! 

XClearArea (display, win, x_des+^, y, 


/ 

if (code->Desorb > code->oldDesorb) 


Light* ='=^->oeJog: code->oldDesorb;_ y _ BoxWidth , heigh t> , 
XFillRectangle (display , win, y _ 


ilse if (code->Desorb < code->oldDesorb) 


height 2 =~code->oldDesorb ^ “^^ D ^°Eoiwidth- 2 , height, 0 ); 

XClearArea (display, win, x_des+/, y, 


if ( code->C02_Ab > code->oldC02_Ab) 


height . 2 ="" code->C 02 ^Ab^- ^>^1^2, y, BoxWidth, height); 

XFillRectangle (display, win, gc, x_ 


glse if (code->C02_Ab < code->oldC02_Ab) 


Light ^Ls“!-“ ! BcxWidth, height, 0); 

XClearArea (display , win, _ 


J 

if ( code->C02_De > code->oldC02_De) 


{ 


LIgh”=Lode->C 02 2 DeL y, BoxWidth , height) 

XFillRectangle (display , win, y _ 


ilse if ( code->C02_De < code->oldC02_De) 


. i q o M rod6“^oldC02 Ds / 
height = J°de->oldC02De - ^e-^ -// BoxWidth, height, o); 
XClearArea (display , wj.h, _ 

/* draw the fan moving */ 
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if 

{ 


} 


(code->sys == 1) 

XClearArea (display, win, *- f “' 

XDrawArc ( display ,win,gc,410, 270, 30, 30, 0,36 

XDrawLine (display , win, gc, 425, 270, , ) » 

XDrawLine (display, win ,gc, 4 10, 285, 440, 285) , 
XDrawLine(display,win,gc,416,275,435,295 , 

XDrawLine (display , win ,gc, 416 , 295, 435, 27b), 


} /* draw_data */ 


C_Height,0) ; 
64) ; 
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